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NOMENCLATURE 


AOA  Angle  of  Attack 

CL  Centerline 

L  Chordwise  length  of  model 

L.E.  Leading  edge 

Mn  Mach  number  normal  to  leading  edge 

U  Velocity  component  in  X  direction 

UINF  Freestream  velocity 

X  Axial  distance  from  nose 

Y  Spanwise  distance  from  centerline 

Z  Distance  normal  to  surface 

a  Angle  of  attack 

On  Angle  of  attack  normal  to  leading  edge 

Alb  leading-edge  sweep  angle 
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1.0  INTRODUCTION 


The  flow  about  modem  high  speed  delta  wing  aircraft  is  characterized  by  the  presence 
of  large  vortices  on  the  leeside  of  the  wing.  In  many  cases,  these  vortices  are  the  primary 
aerodynamic  structure  affecting  the  performance  of  the  aircraft.  The  vortices  form  at  iow  angles 
of  attack,  with  secondary  and  tertiary  vortices  fonning  at  higher  angles  of  attack  due  to  boundary 
layer  separation.  Previous  tests  were  conducted  at  the  Wright  Laboratory  to  obtain  surface 
pressure  profiles,  surface  oil  flow  visualizations,  and  off-body  laser  light  sheet  visualizations 
(Fig.  1).  These  tests  reveal  a  major  change  in  the  supersonic  vortex  structure  above  28°  angle 
of  attack  (AOA)  at  Mach  1 .9.  The  disappearance  of  feeding  sheets,  secondary  vortices,  and 
vortex  shocks  is  observed  in  the  light  she  c  visualizations. 

Laser  velocimetry  (LV)  was  successfully  applied  to  vortex  flows  for  the  subsonic  case  by 
Weissman  (Ref  1)  ,  Yanta  and  Wardlaw  (Ref  2),  Owen  and  Johnson  (Ref  3),  and  Schwind  and 
Mullen  (Ref  4),  among  others.  However,  LV  measurements  have  not  previously  been  reported 
in  the  case  of  a  supersonic  vortex  flowHeld.  Due  to  the  statistical  nature  of  LV  measurements, 
aperiodic,  unsteady  flows  are  biased  by  excessive  turbulence  levels  and  distorted  mean  velocities. 
In  this  test  series,  it  is  critical  to  define  regions  where  the  vortices  burst  so  as  to  preclude 
representing  the  flow  as  highly  turbulent,  but  steady  state. 

Miller  and  Wood  (Ref  5)  present  a  classification  graph  for  flows  over  delta  wings  that 
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Fig.  1.  Laser  Light  Sheet  Visualization  of  Delta  Wing  at  Mach  1.9  and  a 


classify  a  flowfield  into  seven  categories  (Fig.  2):  classical  vortex,  vortex  >vith  shock,  separation 
bubble  with  no  shock,  separation  bubble  with  shock,  shock  with  no  separation,  shock-induced 
separation,  and  no  shock/no  separation.  McMillan,  et  al.  (Ref  6)  verified  the  classification  chart 
using  computational  results  rather  than  the  vapor  screen  flow  visualization  techniques  used  by 
Miller  and  Wood.  The  regions  are  classified  on  a  graph  of  M^  versus  on  where  M^  is  the 
component  of  Mach  number  normal  to  the  leading  edge  and  is  the  angle  of  attack  normal  to 
the  leading  edge  as  given  in  equations  (1)  and  (2): 

Mn  =  M  cos  Ale  (1  +  sin^a  tan^ALe)’^  (1) 

On  =  tan  '(tan  a/cos  Ale)  (2) 

Laser  light  sheet  flow  visualization  of  the  current  75®  delta  wing  shows  that  the  flow  over  this 
model  is  a  classical  vortex  with  a  shock.  The  current  test  extends  the  region  on  the  classification 
graph  representing  "vortex  with  shock"  to  higher  values  of  on- 

For  this  experiment,  a  three-dimensional  LV  system  is  used  to  measure  the  velocity 
components  of  a  supersonic  vortex  flowfield  over  a  simple  sharp-edged  delta  wing  at  Mach  1.9. 
Mean  velocity  vectors  are  presented  for  two  planes  normal  to  the  model  at  three  angles  of  attack. 
Only  a  few  planes  of  data  are  taken,  since  this  study  is  accompanied  by  a  companion 
computational  fluid  dynamic  analysis  of  this  problem  (Ref  7).  The  LV  measurements  are  needed 
for  validation  of  the  CFD  solution,  and  once  the  CFD  solution  is 


Fig.  2.  Classification  Chart  for  Vortices  from  Miller  and  Wood  (1985). 


validated,  the  CFD  solution  will  be  used  as  a  definition  of  the  complete  flowfield.  Angle  of 
attack  effects  are  discussed,  along  with  an  analysis  of  measurement  bias. 
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2.0  TEST  DESCRIPTION 


The  vortical  flow  on  a  sharp-edged  delta  wing  is  measured  at  3  angles  of  attack  (20°,  30°, 
and  35°).  The  leading  edge  sweep  of  the  delta  wing  model  is  75°  and  the  leading  edge  vertex 
angle  is  35°.  A  three-view  sketch  of  the  delta  wing  model  and  a  photograph  of  this  model  are 
shown  in  Figs.  3  and  4.  The  model  length  is  33.91  cm  with  a  thickness  of  1.91  cm.  The  sharp 
edges  have  a  maximum  tip  radius  of  0.025  cm. 

This  vortex  flow  experiment  is  conducted  in  the  Wright  Laboratory’s  Trisonic  Gasdynamics 
Facility  (TGF).  The  TGF  is  a  closed  loop,  continuous  flow  trisonic  wind  tunnel  with  a  0.61  m 
X  0.61  m  (2  ft  X  2  ft)  test  section  for  subsonic  and  supersonic  operations.  For  this  experiment, 
the  TGF  is  operated  at  Mach  1.9  with  a  nominal  stagnation  pressure  of  57.5  kPa  and  stagnation 
temperature  of  310  K.  The  Reynolds  number  for  this  test  is  2.4  x  10*  based  on  the  root  chord 
of  the  model.  The  dynamic  pressure  is  approximately  22.0  kPa. 

A  three-component  laser  velocimeter  (Figs.  5  and  6),  Bragg  shifted  in  the  velocity  component 
normal  to  the  free  steam  and  model  surface,  is  used  to  make  the  measurements.  The  collection 
optics  (Fig.  7)  are  positioned  approximately  10  degrees  off-axis,  limiting  the  effective  loigth  of 
the  probe  volume  to  1.54  mm.  The  flow  is  seeded  with  10-centistoke  silicon  oil  which  is 
introduced  into  the  stagnation  chamber  using  a  specially  designed  seeder.  Out-of-tunnel  particle 
size  measurements  conducted  with  this  seeder  (Ref  8)  show  that  98%  of  the  seed  are  smaller 
than  1.0  ^m  in  diameter,  although  polydisperse. 
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Delta  Wing  Model  in  Trisonic  Gasdynamics  Facility  Test  Section. 


TEST  SECTION 
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Fig.  5.  Three-Component  Laser  Velocimetry  System  for  Delta  Wing  Test. 


Fig.  6.  Transmitting  Optics  for  Laser  Velocimetry  System. 
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The  system  is  configured  in  forward  scatter  mode.  The  transmitting  optics  (Fig.  6)  and 
collection  optics  (Fig.  7)  are  driven  separately  on  two  machinist-style  traverse  tables  with 
accuracy  of  0.0025  mm.  Tunnel  movement  (which  can  be  up  to  1.9  cm  for  a  3-hour  run)  is 
measured  using  a  1-mW  Helium-Neon  laser  and  a  photodiode  position  detecting  system  (Fig. 
8),  with  the  results  sent  to  the  main  computer  and  corrected  coordinates  sent  to  the  traverse 
tables.  The  signals  from  the  5 14-nm-wavelength  green  and  476-nm- wavelength  light  blue  beams 
are  processed  by  a  counter  processor  with  a  1 -nanosecond  clock.  The  signals  from  the  Bragg 
shifted  488  nm  wavelength  blue  beams  are  weaker  and  require  a  Fast  Fourier  Transform  (FFT) 
analysis  to  distinguish  the  signals  from  the  noise;  hence,  the  blue  signals  are  processed  with  a 
hard  wired  processor  based  on  the  FFT  algorithm.  The  counter  processor  and  FFT  analyzer  are 
not  connected;  therefore,  the  multi-component  velocity  data  are  not  restrained  by  a 
predetermined  coincidence  window.  This  allows  for  data  rates  sufficient  to  measure  velocities 
closer  to  the  vortex  core.  The  counterprocessors,  FFT  analyzer,  personal  computer,  and  mini¬ 
computer  are  shown  in  Fig.  9. 

A  schematic  of  the  generic  flowfield  over  the  delta  wing  model  is  shown  in  Fig.  10.  Because 
this  is  a  symmetrical  model,  the  flow  on  only  one-half  of  the  span  is  measured.  The  test  matrices 
consist  of  uniform  grids  on  half  planes  normal  to  the  model  at  40%  and  80%  chord.  Typical 
grid  spacing  is  0.51  cm  by  0.51  cm.  Based  on  previous  laser  light  sheet  and  schlieren  flow 
visualizations,  the  40%  chord  location  is  representative  of  unburst  vortex  flow  and  the  80% 
chord  location  to  represents  burst  vortex  flow  at  the  higher  angles  of  attack. 
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Fig.  8.  Photodiode  Position  Detection  System. 
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Fig.  9.  Data  Collection  System. 
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SCHEMATIC  FLOW  DIAGRAM: 

(1)  SCROLL  VORTEX  SHEET  FORMING  THE  MAIN  VORTEX 

(2)  SECONDARY  VORTEX 

(3)  CENTRAL  ZONE  WITHOUT  VORTEX  FLOW 

(4)  ZONE  INDUCED  BY  THE  MAIN  VORTEX 

(5)  ZONE  INDUCED  BY  THE  SECONDARY  VORTEX 

(6)  ACCUMULATION  ZONE  (AIR  BUBBLE  OR  COATING) 

(7)  REATTACHMENT  LINE 


Fig.  lO.  Schematic  Diagram  of  Flow  about  a  Sharp-Edged  Delta  Wing. 


15 


3.0  MEASUREMENT  ACCURACY 


In  general,  LV  measurements  are  biased  by  many  sources  within  the  optical,  electronic,  data 
processing,  and  seeding  subsystems.  While  some  of  these  sources,  such  as  certain  types  of 
statistical  bias,  may  be  quite  subtle,  other  sources  of  bias  may  dominate  the  measurements.  A 
categorical  breakdown  of  many  of  the  known  sources  of  measurement  error  is  given  by  Jones, 
et  al.  (Ref  9)  and  Edwards  (Ref  10),  which  include  extensive  references  of  LV  bias 
investigations.  For  this  test,  particular  sources  of  measurement  error  appear  to  be  dominant, 
and  help  explain  perturbations  and  voids  within  the  presented  vortex  fields.  By  understanding 
the  sources  and  effects  of  possible  bias,  experimental  error  can  be  distinguished  from  flowfield 
structures,  and  the  data  can  be  used  to  its  full  potential  to  supplement  future  experimental  studies 
or  CFD  validation  efforts. 

One  source  of  bias  is  that  the  LV  measurements  are  time  averaged  to  represent  steady-state 
velocities,  while  the  actual  flowfield  is  not  truly  steady  state.  By  visualizing  cross  sections  of 
the  flow  with  a  laser  light  sheet,  the  vortex  structure  appears  to  be  steady  to  the  naked  eye  or 
to  a  video  camera  with  a  60-frame-per-second  resolution.  However,  by  using  a  high-speed  video 
with  a  1 -microsecond  exposure  time,  the  vortex  contains  structures  which  are  both  periodic  and 
random.  This  inherent  unsteadiness  causes  the  measured  turbulence  intensity  to  be  much  higher 
than  for  a  steady  state  flow.  For  example.  Fig.  1 1  shows  a  representative  turbulence  intensity 
field  with  vertical  and  horizontal  lines  representing  the  magnitudes  of  the  vertical  and  axial 
components  of  turbulence  intensity  with  respect  to  the  freestream  velocity.  Overall,  the 
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Fig,  11,  Turbulence  Intensity  Field  at  x/L  =  40%  for  a  =  30®. 
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turbulence  intensity  in  either  direction  is  a  relatively  uniform  2  to  3  percent  across  the  field, 
except  for  a  few  points  along  the  core  boundary  where  the  axial  values  exceed  20%.  Looking 
at  the  corresponding  histograms  at  these  locations,  shown  in  Fig.  12,  the  flow  at  these  points 
is  bi-modal,  which  indicates  two  distinct  flow  conditions  which  alternate  with  time.  The 
resulting  turbulence  intensity  of  the  combined  flow  is  much  larger  than  the  turbulence  intensities 
of  the  two  individual  modes. 

The  second  important  source  of  bias  within  these  surveys  is  due  to  the  low  signal-to-noise 
ratio  of  the  burst  signal.  Many  factors  contribute  to  the  size  and  clarity  of  the  Doppler  burst, 
including  the  precision  of  the  optical  alignment,  the  transmission  of  light  through  wind  tunnel 
windows,  and  the  quality  and  limitations  of  the  optics  and  electronics  in  general.  During  these 
surveys,  the  data  rate  decreased  substantially  as  the  measurement  location  approached  the  vortex 
core,  as  the  transparency  of  the  windows  is  degraded  by  streaks  of  silicon  oil,  and  as  the  optical 
alignment  drifted  over  a  period  of  several  hours.  Areas  within  the  presented  data  where  changes 
between  neighboring  velocity  vectors  appear  erratic  are  most  likely  effected  by  this  type  of  bias. 

Another  type  of  bias  is  due  to  the  sensitivity  of  the  flowfield  to  slight  changes  in  tunnel 
conditions  during  a  run  and  between  runs.  Most  surveys  required  two  to  three  separate  runs  of 
approximately  5  hours  each.  Despite  efforts  to  precisely  maintain  steady  flow  conditions,  slight 
changes  in  the  flowfield  are  apparent  within  the  presented  velocity  fields.  Rows  or  columns  of 
vectors  with  slightly  different  trends  than  neighboring  rows  and  columns  are  indicative  of  data 
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Fig.  12.  Axial  Velocity  Histograms  around  the  Vortex  Core  Region  at  x/L  =  40%  for  a  =  30“. 


which  are  not  taken  in  sequence  during  a  run,  such  as  certain  rows  of  vectors  within  the  2(f, 
80%  chord  survey.  Areas  where  blocks  of  vectors  within  a  field  exhibit  slightly  different  trends, 
such  as  the  upper  five  rows  of  the  20°,  40%  chord  survey,  represent  data  collected  during 
separate  runs. 

As  LV  measurements  are  applied  to  increasingly  complex,  high  speed  flows,  particle  dynamic 
bias  becomes  a  primary  concern.  Within  the  vortex,  seed  particles  tend  to  be  centrifuged  away 
from  the  core.  This  causes  the  measured  velocity  of  the  seed  to  differ  from  the  fluid  streamlines 
it  is  modeling.  Also,  since  the  larger  particles  are  centrifuged  more  easily  than  smaller  ones, 
the  population  density  and  size  distribution  of  the  seed  vary  radially  through  a  vortex.  A  noted 
decrease  in  data  rate  toward  the  vortex  core  indicates  the  presence  of  fewer,  smaller  particles. 

A  quantitative  analysis  of  the  particle  dynamic  bias  for  the  20°  AOA  flowfield  is  given  by 
Maurice  (Ref  1 1).  By  two  methods  of  analysis,  one  which  models  the  flowfield  by  an  equivalent 
potential  vortex,  and  one  which  uses  a  CFD  flowfield  solution,  the  measurement  error  is  shown 
to  increase  nearly  linearly  with  particle  diameter.  The  actual  particle  size  distribution  of  the 
silicon  oil  is  unknown,  although  out-of-tunnel  measurements  have  shown  98%  of  the  seed  to  be 
less  than  1.0  /im  in  diameter.  The  analytical  analysis  estimates  that  the  seed  must  be  less  than 
0.1  to  0.8  fitn  to  assure  dynamic  bias  of  less  than  3%.  The  analysis  also  shows  that  the  inner 
core  region  did  not  maintain  enough  seed  to  make  measurements,  due  to  the  inability  of  particles 
to  remain  in  the  core  region  as  the  vortex  is  formed,  and  due  to  particles  outside  of  the  viscous 
inner  core  being  centrifuged  as  they  travel  downstream.  Consequently,  a  specific  estimation  of 
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the  particle  dynamic  bias  requires  precise  knowledge  of  the  particle  diameter,  and  maintaining 
seed  within  the  inner  core  cannot  be  accomplished  by  injecting  the  particles  into  the  freestream 
flow. 


Attempts  to  seed  through  the  model  into  the  vortex  at  its  creation  result  in  significant  changes 
in  the  vortex  structure.  As  a  result,  this  is  no  longer  a  nonintrusive  measurement  and  could  not 
be  correlated  with  computational  analyses;  hence,  no  data  are  presented  from  this  seeding 
technique.  Future  tests  will  be  conducted  with  tailored  flow  seeding  through  the  model. 
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4.0  RESULTS 


In  general,  1500  samples  are  taken  at  each  test  location.  However,  at  some  locations, 
especially  near  the  virtually  unseeded  vortex  core,  data  rates  are  too  slow  to  acquire  the  1500 
samples  and  the  data  collection  is  terminated  after  3  minutes.  Model  interference  with  the  laser 
beams  limits  access  to  some  test  locations  near  the  surface  of  the  model.  These  difficulties 
complicate  the  problem  of  defining  the  secondary  vortices. 

Two-dimensional  velocity  fields  at  20°,  30°,  and  35°  AOA  for  40%  and  80%  chord  stations 
are  shown  in  Figs.  13-17.  The  small  isolated  holes  in  the  flow  field  represent 
instrumentation/computer  errors.  The  large  empty  areas  near  the  model  surface  are  areas  in, 
or  around  the  vortex  core  where  the  LV  seed  material  is  centrifuged  away  from  the  vortex  core. 
Also,  there  is  a  region  along  the  surface  where  the  laser  beams  intersect  the  model  prior  to  the 
formation  of  the  probe  volume,  restricting  LV  measurements  in  the  primary  vortex.  Because 
of  these  two  limitations,  the  secondary  and  tertiary  vortices  cannot  be  measured  without  set-up 
modifications  that  are  being  considered  for  further  studies. 

At  20°  AOA,  previous  experiments  and  laser  light  sheet  flow  visualizations  indicate  that  there 
is  not  any  bursting  of  the  vortices  (Ref  12).  Hence,  the  LV  measurements  provide  more 
accurate  statistical  mean  velocity  fields  at  20“  AOA  than  at  30°  or  35°  AOA.  At  40%  chord  the 
downward  flow  along  the  centerline  indicates  flow  reattachment  as  illustrated  in  the  schematic 
flow  diagram  (Fig.  10).  The  vertical  vectors  along  the  centerline  in  Figs.  13  and  14  clearly 
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Fig.  13.  Two-Dimensional  Velocity  Field  at  x/L  =  40%  for  a  =  20°. 
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Fig.  14.  Two-Dimensional  Velocity  Field  at  x/L  =  80%  for  a  =  20®. 
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Fig.  15.  Two-Dimensional  Velocity  Field  at  x/L  =  40%  for  a  -  30®. 
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Fig.  16.  Two-Dimensional  Velocity  Field  at  x/L  =  80%  for  o  =  30®. 
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Fig.  17.  Two-Dimensional  Velocity  Field  at  x/L  =  40%  for  a  =  35". 
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show  the  symmetry  of  the  flow.  The  quasi-randomness  of  the  vectorfield  at  35®  AOA,  40% 
chord  represents  the  vortex  bursting.  At  80%  chord,  the  vortex  bursting  is  too  strong  for  this 
LV  system  to  measure. 


A  representative  velocity  histogram  of  a  single  measurement  component  at  a  single  location 
is  shown  in  Fig.  18.  After  1500  individual  point  measurements  are  made,  the  data  are  edited 
by  eliminating  the  points  which  fall  beyond  3  standard  deviations  of  the  mean.  At  this  location 
30  points  are  eliminated  by  this  form  of  editing.  The  remaining  1470  points  yield  a  mean 
velocity  of  483.2  m/s  with  a  standard  deviation  of  22.24  m/s.  The  ideal  Gaussian  distribution 
curve  with  the  same  mean  and  variance  is  superimposed  over  the  data.  At  this  point  the 
skewness  coefficient  is  -0.39,  relative  to  a  coefficient  of  zero  for  the  ideal  Gaussian  curve.  The 
kurtosis  coefficient  at  this  point  is  2.82,  compared  to  an  ideal  value  of  3.00. 

A  2-D  velocity  field  for  the  20°  AOA,  80%  chord  survey  where  the  data  are  rotated  to  a 
coordinate  system  which  is  axial  and  normal  to  the  trajectory  of  the  vortex  core  is  presented  in 
Fig.  19.  Since  the  model  is  at  a  20  degree  angle-of-attack  relative  to  the  tunnel,  and  since  the 
vortex  core  grows  as  it  moves  downstream  along  the  model,  the  origin  of  the  vortex  is  trajected 
at  a  lesser  angle-of-attack  than  the  wing.  By  extending  "tails"  from  the  smooth,  unbiased 
velocity  vectors  as  shown  in  the  figure,  a  point  of  convergence  defines  the  vortex  origin.  For 
a  transformation  angle  greater  than  8®  (towards  the  20®  AOA),  the  lines  converge  to  create  a 
surface  concave  to  the  right.  For  an  angle  of  less  than  8®,  the  surface  is  concave  towards  the 
left.  Since  a  small  perturbation  in  vector  orientation  creates  a  large  change  in  the 
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Fig.  18.  Typical  Velocity  Histogram  of  Axial  Component  of  Velocity. 


a  =  20  “  x/L  =  0.8 
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Fig.  19.  Transformed  Two-Dimensional  Velocity  Field  Estimating  Lxx^tion  of  Vortex 
Origin  at  x/L  =  80%  for  a  =  20". 
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position  of  its  tail,  this  method  requires  very  smooth,  unbiased  data  to  be  effective.  Slight 
changes  in  flowfleld  conditions  over  a  long  tunnel  run  are  evidenced  by  differ»it  groups  of 
vectors  converging  to  different  origins  within  the  same  velocity  field.  Knowing  the  position  of 
the  vortex  origin  allows  each  measurement  location  to  be  related  to  an  ideal  potential  vortex 
strength.  For  this  case,  the  potential  vortex  circulation  increases  from  approximately  25  mVs 
at  the  centerline  to  100  mVs  at  the  outer  edge  of  the  measurement  field. 

A  three-dimensional  representation  of  the  axial  velocity  component  on  the  40%  chord  plane 
at  20®  AOA  (Fig.  20)  shows  the  momentum  loss  in  the  vortex  region  of  the  flowfield  as 
evidenced  by  the  decreasing  magnitudes  of  the  axial  velocity.  Figs.  21  and  22  quantitatively 
show  the  velocity  trends  in  this  flowfield. 

U/UINF  versus  Y/L  for  2JL  =  0.(X)7  to  0.075  is  shown  in  Fig.  21.  In  general,  the  ratio  of 
U/UINF  increases  as  the  distance  normal  to  the  model  surface  increases  at  all  locations  along 
the  span.  At  the  centerline  of  the  model,  the  normalized  axial  velocity  increases  from  near  0.9 
to  approximately  1.13  as  Z/L  increases  from  0.007  to  0.052,  and  then  the  normalized  axial 
velocity  decreases  to  1  .(X)  as  Z/L  increases  to  0. 142  (the  far  field).  Considering  the  ftrst  full 
set  of  data  at  Z/L  =  0.022,  it  can  be  seen  that  the  velocity  is  highest  at  the  centerline  and 
between  87%  and  93%  span  and  that  this  last  location  group  represents  a  stagnation  line  between 
the  primary  and  secondary  vortices  on  an  oil  flow  simulation.  This  compares  qualitatively  to 
the  stagnation  line  on  the  simulated  oil  traces  of  Webster  and  Shang  (Ref  13).  At  approximately 
50%  span,  near  the  spanwise  location  of  the  vortex  core,  the  normalized  velocities  increase  from 
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Fig,  20.  Three-Dimensional  Representation  of  Axial  Velocity  Component  on  the  x/L  =  40% 
Plane  for  a  =  20“. 
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Fig.  21.  Normalized  Axial  Velocity  Component  vs.  Span-wise  Distance  on  Model  for 
z  =  0.25  cm  to  z  =  2.54  cm. 
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Fig.  22.  Normalized  Axial  Velocity  Component  vs.  Span-wise  Distance  on  Model  for 
z  =  2.79  cm  to  z  4.83  cm. 


0.44  to  0.98  as  Z/L  increases  from  0.008  to  0.07S. 

U/UINF  versus  Y/L  for  Z/L  =  0.075  to  0. 142  is  shown  in  Fig.  22.  The  band  of  normalized 
velocity  ranges  from  1.0  to  1.1  at  the  coiterline  of  the  model.  The  shaded  band  shows  the 
velocity  trends  outside  of  the  vortex  core.  This  compares  qualitatively  with  the  results  of  Gad* 
el-Hak  and  Blackwelder  (Ref  14,  IS)  for  the  subsonic  delta  wing  flowfield  near  the  model 
surface.  The  velocity  decreases  as  the  spanwise  distance  from  the  centerline  increases,  that  rises 
significantly  at  the  leading  edge,  and  finally  tapers  off  to  approximately  95%  of  freestieam 
velocity  just  off  the  edge  of  the  model. 
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5.0  CONCLUSIONS 


An  experimental  study  was  conducted  to  measure  velocities  in  the  vortices  generated  by 
a  sharp  leading  edge  75°  delta  wing  at  Mach  1.9  for  angles  of  attack  up  to  35°.  Three- 
dimrasional  laser  velocimeter  measurements  were  made  to  define  the  vortex  structure  prior  to 
bursting  at  several  chord  locations  as  part  of  an  ongoing  and  developing  project  for 
computational  fluid  dynamics  code  verification  and  validation.  Qualitative  comparisons  between 
experimental  data  and  CFD  results  show  good  agreement  on  defining  the  region  between  the 
primary  and  secondary  vortices.  However,  the  CFD  solutions  revealed  no  evidence  of  vortex 
breakdown  which  was  visualized  using  a  laser  light  sheet  and  seen  at  30°  angle  of  attack  and 
80%  chord.  An  analysis  of  the  light  sheet  results  is  ongoing.  Comparison  of  test  characteristics 
and  results  with  the  vortex  classification  plot  of  Miller  and  Wood  (Ref  5)  show  that  the 
classification  graph  can  be  extended  to  higher  normalized  angle  of  attack  in  the  classical  vortex 
region.  Measurements  of  the  vortex  core  are  shown  to  be  limited  by  the  lack  of  seed  material 
entrained  in  the  vortex  and  by  the  centrifugal  effect  of  the  vortex  on  the  seed  material. 
Additional  experiments  are  planned  with  tailored  flow  seeding  through  the  model  and  enhanced 
electronic  equipment  to  further  define  the  vortex  structure,  including  secondary  and  tertiary 
vortices,  turbulence,  and  stability. 
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6.0  RECOMMENDATIONS 


Recommendations  for  this  program  include  improving  seeding  techniques  and  a  better 
mixing  of  electronic  equipment.  As  a  result  of  this  experimratal  study,  the  tedmiques  of 
particle  dynamics  was  developed  by  Maurice  (Ref  11,  see  Appendix  B).  With  this  additional 
analysis  technique,  experimental  studies  with  different  seeding  techniques  and  seed  particles  can 
be  conducted  to  optimize  particle  tracking  and  bias.  Also,  the  currrat  experiments  showed  that 
use  of  a  Fast  Fourier  Transform  processor  for  the  weaker  laser  velocimetry  components  allowed 
those  measurements  which  were  unobtainable  widi  the  counterprocessors.  At  the  time  of  these 
experiments,  only  one  FFT  processor  could  reasonably  be  used,  but  now  the  equipment  exists 
to  perform  this  study  using  three  FFT  processors  to  collect  coincident  data.  With  coincident 
data,  additional  turbulence  quantities  can  be  determined  from  the  data.  Also,  much  more 
thorough  data  collection  could  be  made  if  unlimited  wind  tunnel  time  was  available  (but  that’s 
just  an  experimenter  dreaming). 
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APPENDIX  A:  QUANTITATIVE  RESULTS 


This  appendix  contains  the  numerical  results  of  the  laser  velocimetry  measurements  for 
20"  AOA,  40%  chord;  20"  AOA,  80%  chord;  30"AOA,  40%  chord;  30"  AOA,  80%  chord;  and 
35",  40%  chord.  The  coordinates  and  velocities  are  presented  in  the  model  coordinate  system 
as  shown  in  Fig.A-1.  The  coordinates  x,  y,  and  z  are  given  in  meters  and  the  respective 
velocity  components  u,  v,  and  w  are  given  in  meters  per  second.  Any  locations  where  no  results 
are  given  were  locations  where  the  data  are  invalid,  not  acquired  due  to  electronic  errors,  or  not 
acquired  due  to  particle  dynamics. 
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Table  A-1:  Experimental  Data  for  20“  Angle  of  Attack,  40%  Chord. 


X 

y 

z 

u 

V 

w 

(m) 

(m) 

(m) 

(m/s) 

(m/s) 

(m/s) 

-0.1357 

0.0508 

0.0025 

137.75 

0.19 

52.25 

-0.1357 

0.0508 

0.0076 

137.59 

-0.48 

52.97 

-0.1357 

0.0508 

0.0127 

138.42 

-3.80 

54.02 

-0.1357 

0.0508 

0.0178 

139.34 

-4.75 

54.98 

-0.1357 

0.0508 

0.0229 

140.44 

-6.44 

55.86 

-0.1357 

0.0508 

0.0279 

139.66 

-9.89 

55.70 

-0.1357 

0.0508 

0.033 

139.62 

-4.02 

55.62 

-0.1357 

0.0508 

0.0381 

142.26 

-6.71 

56.45 

-0.1357 

0.0508 

0.0432 

144.66 

-14.09 

57.40 

-0.1357 

0.0508 

0.0483 

143.18 

-11.55 

56.84 

-0.1357 

0.0508 

0.0508 

143.36 

-13.01 

56.99 

-0.1357 

0.0457 

0.0025 

136.47 

-2.02 

51.65 

-0.1357 

0.0457 

0.0076 

125.50 

12.93 

48.36 

-0.1357 

0.0457 

0.0127 

130.26 

5.03 

51.80 

-0.1357 

0.0457 

0.0178 

138.50 

-4.40 

55.93 

-0.1357 

0.0457 

0.0229 

140.11 

-6.10 

56.86 

-0.1357 

0.0457 

0.0279 

140.15 

-7.30 

56.77 

-0.1357 

0.0457 

0.033 

139.59 

-5.85 

56.14 

-0.1357 

0.0457 

0.0381 

142.79 

-7.81 

57.14 

-0.1357 

0.0457 

0.0432 

145.71 

-15.66 

58.17 

-0.1357 

0.0457 

0.0483 

145.14 

-14.68 

57.82 

-0.1357 

0.0457 

0.0508 

146.25 

-19.04 

58.17 

-0.1357 

0.0432 

0.0025 

108.31 

45.02 

33.63 

-0.1357 

0.0432 

0.0051 

114.26 

44.66 

41.72 

-0.1357 

0.0432 

0.0076 

118.52 

39.66 

48.65 

43 


nn 

y 

z 

u 

V 

w  1 

(m) 

(m) 

(m/s) 

(m/s) 

(m/s)  1 

-0.1357 

0.0432 

0.0102 

121.49 

32.16 

53.55  1 

-0.1357 

0.0432 

0.0127 

123.61 

24.82 

57.08  1 

-0.1357 

0.0432 

0.0152 

125.15 

17.95 

57.54  1 

-0.1357 

0.0432 

0.0178 

126.35 

12.43 

57.77  1 

-0.1357 

0.0432 

0.0203 

128.52 

7.82 

58.57  1 

-0.1357 

0.0432 

0.0229 

133.71 

1.88 

58.05 

-0.1357 

0.0432 

0.0254 

151.52 

-14.79 

63.17 

-0.1357 

0.0432 

0.0279 

144.95 

-13.55 

60.78 

-0.1357 

0.0432 

0.0305 

143.67 

-12.82 

60.31 

-0.1357 

0.0432 

0.033 

145.45 

-13.97 

60.33 

-0.1357 

0.0432 

0.0356 

151.90 

-16.51 

62.98 

-0.1357 

0.0432 

0.0381 

153.95 

-15.87 

62.58 

-0.1357 

0.0432 

0.0407 

152.64 

-20.42 

61.85 

-0.1357 

0.0432 

0.0432 

153.80 

-20.22 

61.99 

-0.1357 

0.0407 

0.0025 

119.34 

19.96 

42.69 

-0.1357 

0.0407 

0.0051 

118.10 

39.48 

59.92 

-0.1357 

0.0407 

0.0076 

123.21 

22.03 

63.73 

-0.1357 

0.0407 

0.0102 

126.62 

-4.44 

65.84 

-0.1357 

0.0407 

0.0127 

127.48 

-19.70 

77.66 

-0.1357 

0.0407 

0.0152 

131.66 

-21.06 

81.82 

-0.1357 

0.0407 

0.0178 

132.63 

-13.25 

79.34 

-0.1357 

0.0407 

0.0203 

135.77 

-15.60 

74.93 

-0.1357 

0.0407 

0.0229 

137.52 

-17.66 

71.67 

-0.1357 

0.0407 

0.0254 

138.38 

-16.29 

70.13 

-0.1357 

0.0407 

0.0279 

141.44 

-16.09 

69.11 

-0.1357 

0.0407 

0.0305 

143.47 

-14.78 

66.44 

44 


X 

y 

z 

u 

V 

w 

(m) 

(m) 

(tn) 

(m/s) 

(m/s) 

(m/s) 

-0.1357 

0.0407 

0.033 

145.05 

-15.48 

65.30 

-0.1357 

0.0407 

0.0356 

150.51 

-27.23 

66.89 

-0.1357 

0.0407 

0.0381 

147.75 

-17.96 

63.71 

-0.1357 

0.0407 

0.0407 

149.24 

-18.00 

63.14 

-0.1357 

0.0407 

0.0432 

147.40 

-18.55 

62.05 

-0.1357 

0.0407 

0.0457 

148.63 

-16.30 

61.96 

-0.1357 

0.0407 

0.0483 

149.53 

-19.37 

62.19 

-0.1357 

0.0407 

0.0508 

152.05 

-26.72 

62.50 

-0.1357 

0.0381 

0.0025 

99.49 

35.96 

39.11 

-0.1357 

0.0381 

0.0051 

117.01 

38.12 

48.95 

^.1357 

0.0381 

0.0076 

139.57 

3.04 

58.84 

-0.1357 

0.0381 

0.0102 

140.43 

-2.02 

61.57 

-0.1357 

0.0381 

0.0127 

142.08 

-16.29 

63.20 

-0.1357 

0.0381 

0.0152 

143.52 

-33.48 

63.20 

-0.1357 

0.0381 

0.0178 

139.91 

-27.31 

62.98 

-0.1357 

0.0381 

0.0203 

141.75 

-32.05 

71.89 

-0.1357 

0.0381 

0.0229 

139.80 

-22.69 

76.14 

-0.1357 

0.0381 

0.0254 

140.34 

-24.55 

73.22 

-0.1357 

0.0381 

0.0279 

141.19 

-16.08 

70.32 

-0.1357 

0.0381 

0.0305 

143.63 

-16.32 

69.33 

-0.1357 

0.0381 

0.033 

144.91 

-18.21 

66.40 

-0.1357 

0.0381 

0.0356 

146.96 

-8.28 

61.58 

-0.1357 

0.0381 

0.0381 

147.73 

-11.47 

63.44 

-0.1357 

0.0381 

0.0407 

148.49 

-12.74 

63.11 

-0.1357 

0.0381 

0.0432 

148.98 

-13.62 

62.44 

-0.1357 

0.0381 

0.0457 

149.19 

-12.44 

62.07 
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X 

y 

z 

u 

V 

w  1 

(m) 

(m) 

(m) 

(m/s) 

(m/s) 

(m/s)  1 

1  -0.1357 

0.0381 

0.0483 

149.13 

-12.01 

61.53  1 

1  -0.1357 

0.0381 

0.0508 

148.52 

-10.06 

61.00 

1  -0.1357 

0.0356 

0.0051 

115.91 

13.53 

44.94 

1  -0.1357 

0.0356 

0.0076 

137.88 

-13.44 

55.93 

-0.1357 

0.0356 

0.0102 

140.27 

-2.95 

61.22 

-0.1357 

0.0356 

0.0127 

141.80 

-1.25 

64.47 

-0.1357 

0.0356 

0.0152 

143.17 

-3.49 

66.34 

1  -0.1357 

0.0356 

0.0178 

143.45 

-31.39 

67.60 

-0.1357 

0.0356 

0.0203 

141.40 

-32.66 

68.99 

-0.1357 

0.0356 

0.0229 

140.40 

-28.31 

73.16 

-0.1357 

0.0356 

0.0254 

139.89 

-23.58 

76.41 

-0.1357 

0.0356 

0.0279 

140.20 

-22.04 

72.69 

-0.1357 

0.0356 

0.0305 

140.55 

-8.09 

70.70 

-0.1357 

0.0356 

0.033 

142.82 

-8.30 

67.02 

-0.1357 

0.0356 

0.0356 

145.33 

-9.14 

64.86 

-0.1357 

0.0356 

0.0381 

148.56 

-20.04 

64.36 

-0.1357 

0.0356 

0.0407 

150.97 

-25.05 

65.30 

-0.1357 

0.0356 

0.0432 

151.06 

-26.00 

63.91 

-0.1357 

0.0356 

0.0457 

150.16 

-18.59 

62.27 

-0.1357 

0.0356 

0.0483 

149.41 

-12.01 

61.37 

-0.1357 

0.0356 

0.0508 

149.46 

-12.00 

61.16 

-0.1357 

0.033 

0.0076 

128.76 

7.27 

49.34 

-0.1357 

0.033 

0.0102 

131.97 

-10.23 

52.86 

-0.1357 

0.033 

0.0127 

136.80 

-8.10 

59.35 

-0.1357 

0.033 

0.0152 

128.32 

16.30 

59.85 

-0.1357 

0.033 

0.0178 

144.37 

-26.94 

67.09 

46 


X 

(m) 

y 

(m) 

z 

(m) 

u 

(m/s) 

V 

(m/s) 

w 

(m/s) 

-0.1357 

0.033 

0.0203 

147.73 

-52.75 

70.49 

-0.1357 

0.033 

0.0229 

150.17 

-57.29 

74.37 

-0.1357 

0.033 

0.0254 

148.52 

-50.43 

76.09 

-0.1357 

0.033 

0.0279 

146.95 

-38.08 

75.97 

-0.1357 

0.033 

0.0305 

146.75 

-38.08 

74.17 

-0.1357 

0.033 

0.033 

146.88 

-27.69 

70.13 

-0.1357 

0.033 

0.0356 

149.00 

-28.71 

68.33 

-0.1357 

0.033 

0.0381 

146.27 

-22.59 

64.36 

-0.1357 

0.033 

0.0407 

152.15 

-30.96 

61.88 

-0.1357 

0.033 

0.0432 

150.96 

-30.20 

61.13 

-0.1357 

0.033 

0.0457 

156.47 

-23.28 

61.93 

-0.1357 

0.033 

0.0483 

153.43 

-33.03 

62.61 

-0.1357 

0.033 

0.0508 

153.94 

-32.03 

62.02 

-0.1357 

0.0305 

0.0076 

124.55 

22.83 

32.77 

-0.1357 

0.0305 

0.0102 

133.61 

11.02 

36.72 

-0.1357 

0.0305 

0.0127 

135.41 

0.60 

43.41 

-0.1357 

0.0305 

0.0152 

144.19 

-9.25 

57.13 

-0.1357 

0.0305 

0.0178 

146.08 

-14.33 

59.84 

-0.1357 

0.0305 

0.0203 

142.93 

-6.13 

64.99 

-0.1357 

0.0305 

0.0229 

156.01 

-39.50 

71.20 

-0.1357 

0.0305 

0.0254 

154.59 

-37.78 

74.01 

-0.1357 

0.0305 

0.0305 

153.49 

-27.45 

67.13 

-0.1357 

0.0305 

0.033 

136.11 

7.13 

56.61 

-0.1357 

0.0305 

0.0356 

145.2/ 

-30.62 

68.53 

-0.1357 

0.0305 

0.0381 

146.24 

-21.22 

65.04 

-0.1357 

0.0305 

0.0407 

150.18 

-24.68 

60.90 
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X 

y 

z 

u 

V 

w 

(m) 

(m) 

(m) 

(m/s) 

(m/s) 

(m/s) 

-0.1357 

0.0305 

0.0432 

151.33 

-26.21 

61.21 

-0.1357 

0.0305 

0.0457 

152.03 

-28.37 

61.22 

-0.1357 

0.0305 

0.0483 

152.39 

-27.72 

61.29 

-0.1357 

0.0305 

0.0508 

153.70 

-32.59 

61.72 

-0.1357 

0.0279 

0.0076 

95.05 

44.02 

24.02 

-0.1357 

0.0279 

0.0102 

119.68 

17.72 

29.18 

-0.1357 

0.0279 

0.0127 

127.39 

5.45 

33.79 

-0.1357 

0.0279 

0.0152 

134.93 

-18.15 

52.66 

-0.1357 

0.0279 

0.0178 

135.78 

-22.17 

54.39 

-0.1357 

0.0279 

0.0203 

140.52 

-16.96 

61.55 

-0.1357 

0.0279 

0.0229 

143.49 

-41.37 

64.87 

-0.1357 

0.0279 

0.0254 

147.65 

-52.98 

66.12 

-0.1357 

0.0279 

0.0279 

150.47 

-42.94 

66.23 

-0.1357 

0.0279 

0.0305 

153.51 

-49.91 

66.47 

-0.1357 

0.0279 

0.033 

152.41 

-47.95 

65.00 

-0.1357 

0.0279 

0.0356 

149.16 

-27.40 

64.74 

-0.1357 

0.0279 

0.0381 

148.49 

-16.79 

63.33 

-0.1357 

0.0279 

0.0407 

150.57 

-28.79 

61.18 

-0.1357 

0.0279 

0.0432 

150.68 

-28.80 

61.14 

-0.1357 

0.0279 

0.0457 

152.06 

-28.55 

61.35 

-0.1357 

0.0279 

0.0483 

152.39 

-27.07 

61.23 

-0.1357 

0.0279 

0.0508 

153.23 

-33.18 

61.43 

-0.1357 

0.0254 

0.0127 

108.28 

33.73 

33.28 

-0.1357 

1 

0.0254 

0.0152 

126.09 

-2.66 

41.23 

-0.1357 

0.0254 

0.0178 

128.80 

-14.02 

45.61 

-0.1357 

0.0254 

0.0203 

137.72 

-27.29 

55.03 
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X 

(m) 

y 

(m) 

z 

(m) 

u 

(m/s) 

V 

(m/s) 

w 

(m/s) 

-0.1357 

0.0254 

0.0229 

148.34 

-30.13 

62.65 

-0.1357 

0.0254 

0.0254 

159.63 

-56.96 

75.79 

-0.1357 

0.0254 

0.0279 

158.04 

-49.90 

76.00 

-0.1357 

0.0254 

0.0305 

156.91 

-39.88 

73.92 

-0.1357 

0.0254 

0.033 

155.45 

-29.78 

70.89 

-0.1357 

0.0254 

0.0356 

154.05 

-21.75 

68.15 

-0.1357 

0.0254 

0.0381 

148.80 

-7.06 

64.23 

-0.1357 

0.0254 

0.0407 

155.41 

-30,55 

62.66 

-0.1357 

0.0254 

0.0432 

156.41 

-30.18 

61.98 

-0.1357 

0.0254 

0.0457 

156.43 

-23.47 

61.19 

-0.1357 

0.0254 

0.0483 

155.98 

-24.12 

60.55 

-0.1357 

0.0254 

0.0508 

156.05 

-20.06 

59.97 

-0.1357 

0.0229 

0.0076 

86.20 

34.77 

19.91 

-0.1357 

0.0229 

0.0102 

108.94 

34.26 

18.78 

-0.1357 

0.0229 

0.0127 

106.75 

26.68 

19.77 

-0.1357 

0.0229 

0.0152 

118.00 

9.97 

34.31 

-0.1357 

0.0229 

0.0178 

145.49 

-40.11 

47.77 

-0.1357 

0.0229 

0.0203 

136.51 

-24.61 

51.45 

-0.1357 

0.0229 

0.0229 

146.52 

-43.38 

61.86 

-0.1357 

0.0229 

0.0254 

157.93 

-51.69 

65.39 

-0.1357 

0.0229 

0.0279 

160.52 

-55.41 

71.75 

-0.1357 

0.0229 

0.0305 

158.81 

-43.58 

70.27 

-0.1357 

0.0229 

0.033 

157.33 

-33.36 

67.98 

-0.1357 

0.0229 

0.0356 

156.10 

-26.08 

66.55 

-0.1357 

0.0229 

0.0381 

151.41 

-22.25 

62.35 

-0.1357 

0.0229 

0.0407 

155.59 

-32.72 

62.16 

X 

y 

z 

u 

V 

w 

(m) 

(m) 

(m) 

(m/s) 

(m/s) 

(m/s) 

-0.1357 

0.0229 

0.0432 

156.57 

-32.24 

61.38 

-0.1357 

0.0229 

0.0457 

156.69 

-28.35 

60.72 

-0.1357 

0.0229 

0.0483 

156.37 

-22.14 

60.34 

-0.1357 

0.0229 

0.0508 

156.15 

-23.48 

59.76  1 

-0.1357 

0.0203 

0.0076 

66.92 

6.80 

5.40  1 

-0.1357 

0.0203 

0.0102 

92.62 

30.39 

13.92  1 

-0.1357 

0.0203 

0.0127 

95.51 

16.00 

17.39  1 

-0.1357 

0.0203 

0.0152 

102.27 

13.05 

20.75  1 

1  -0.1357 

0.0203 

0.0178 

123.79 

-2.15 

35.74 

-0.1357 

0.0203 

0.0203 

130.27 

-25.28 

44.40 

-0.1357 

0.0203 

0.0229 

139.34 

-35.83 

54.21 

-0.1357 

0.0203 

0.0254 

152.45 

-39.85 

59.89 

-0.1357 

0.0203 

0.0279 

162.30 

-55.66 

65.76 

-0.1357 

0.0203 

0.0305 

160.41 

-43.58 

64.72 

-0.1357 

0.0203 

0.033 

159.05 

-35.10 

63.75 

-0.1357 

0.0203 

0.0356 

156.73 

-26.48 

63.51 

-0.1357 

0.0203 

0.0381 

154.45 

-18.89 

62.90  1 

-0.1357 

0.0203 

0.0407 

156.01 

-30.83 

61.38 

-0.1357 

0.0203 

0.0432 

155.90 

-29.80 

60.63 

-0.1357 

0.0203 

0.0457 

156.54 

-22.87 

60.16 

-0.1357 

0.0203 

0.0483 

156.34 

-24.64 

59.76 

-0.1357 

0.0203 

0.0508 

156.34 

-22.97 

59.30 

-0.1357 

0.0178 

0.0127 

87.36 

21.48 

12.51 

-0.1357 

0.0178 

0.0254 

149.66 

-32.93 

58.33 

-0.1357 

0.0178 

0.0279 

162.95 

-49.64 

59.54 

-0.1357 

0.0178 

0.0305 

161.76 

-39.50 

57.22 
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X 

y 

z 

u 

V 

w 

(m) 

(m) 

(m) 

(m/s) 

(m/s) 

(m/s) 

-0.1357 

0.0178 

0.033 

160.05 

-33.67 

58.33 

-0.1357 

0.0178 

0.0356 

157.87 

-26.36 

59.31 

-0.1357 

0.0178 

0.0381 

155.70 

-20.25 

60.11 

-0.1357 

0.0178 

0.0407 

156.51 

-32.69 

60.58 

-0.1357 

0.0178 

0.0432 

156.16 

-29.55 

59.86 

-0.1357 

0.0178 

0.0457 

156.63 

-28.63 

59.55 

-0.1357 

0.0178 

0.0483 

156.21 

-25.65 

59.17 

-0.1357 

0.0178 

0.0508 

156.02 

-26.38 

58.53 

-0.1357 

0.0152 

0.0279 

166.49 

-39.64 

41.94 

-0.1357 

0.0152 

0.0305 

163.94 

-33.25 

47.60 

-0.1357 

0.0152 

0.033 

161.33 

-28.99 

53.01 

-0.1357 

0.0152 

0.0356 

159.71 

-25.37 

55.85 

-0.1357 

0.0152 

0.0381 

156.52 

-27.82 

57.48 

-0.1357 

0.0152 

0.0407 

156.87 

-32.92 

59.38 

-0.1357 

0.0152 

0.0432 

155.69 

-28.86 

58.81 

-0.1357 

0.0152 

0.0457 

156.09 

-27.58 

58.66 

-0.1357 

0.0152 

0.0483 

155.88 

-27.40 

58.31 

-0.1357 

0.0152 

0.0508 

156.08 

-23.44 

58.06 

-0.1357 

0.0127 

0.0178 

112.85 

2.22 

20.67 

-0.1357 

0.0127 

0.0305 

166.72 

-25.75 

34.12 

-0.1357 

0.0127 

0.033 

163.16 

-22.38 

44.13 

-0.1357 

0.0127 

0.0356 

160.18 

-20.16 

50.80 

-0.1357 

0.0127 

0.0381 

158.30 

-17.20 

53.80 

-0.1357 

0.0127 

0.0407 

156.88 

-29.22 

58.09 

-0.1357 

0.0127 

0.0432 

156.77 

-28.04 

57.87 

-0.1357 

0.0127 

0.0457 

156.54 

-26.86 

57.97 
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X 

y 

z 

u 

V 

w 

(m) 

(m) 

(m) 

(m/s) 

(m/s) 

(m/s) 

-0.1357 

0.0127 

0.0483 

156.62 

-25.93 

57.83 

-0.1357 

0.0127 

0.0508 

156.31 

-23.21 

57.47 

-0.1357 

0.0102 

0.0254 

160.46 

-19.10 

42.08 

-0.1357 

0.0102 

0.0279 

169.21 

-18.81 

21.36 

-0.1357 

0.0102 

0.0305 

166.91 

-17.78 

28.75 

-0.1357 

0.0102 

0.033 

164.85 

-15.46 

35.11 

-0.1357 

0.0102 

0.0356 

161.57 

-15.18 

44.79 

-0.1357 

0.0102 

0.0381 

158.51 

-13.22 

51.29 

-0.1357 

0.0102 

0.0407 

156.64 

-29.69 

57.19 

-0.1357 

0.0102 

0.0432 

156.53 

-28.47 

56.80 

-0.1357 

0.0102 

0.0457 

156.78 

-27.49 

57.24 

-0.1357 

0.0102 

0.0508 

156.46 

-23.24 

56.83 

-0.1357 

0.0076 

0.0254 

170.58 

-10.42 

10.36 

-0.1357 

0.0076 

0.0279 

169.09 

-12.02 

18.34 

-0.1357 

0.0076 

0.0305 

167.49 

-12.64 

25.19 

-0.1357 

0.0076 

0.033 

165.95 

-11.52 

31.30 

-0.1357 

0.0076 

0.0356 

164.08 

-10.82 

36.73 

-0.1357 

0.0076 

0.0381 

159.37 

-8.50 

48.16 

-0.1357 

0.0076 

0.0407 

156.91 

-31.03 

56.02 

-0.1357 

0.0076 

0.0432 

156.82 

-27.03 

56.18 

-0.1357 

0.0076 

0.0457 

156.89 

-24.88 

56.22 

-0.1357 

0.0076 

0.0483 

156.73 

-21.36 

56.43 

-0.1357 

0.0076 

0.0508 

156.64 

-20.24 

56.13 

-0.1357 

0.0051 

0.0076 

109.74 

-4.06 

-15.36 

-0.1357 

0.0051 

0.0102 

109.46 

-9.52 

-12.78 

-0.1357 

0.0051 

0.0127 

146.62 

-33.11 

4.18 
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X 

(m) 

r  y 

(m) 

z 

(m) 

u 

(m/s) 

V 

(m/s) 

w 

(m/s) 

-0.1357 

0.0051 

0.0152 

112.68 

-14.94 

-1.96 

-0.1357 

0.0051 

0.0178 

114.84 

-4.57 

3.04 

-0.1357 

0.0051 

0.0203 

173.00 

-1.43 

-9.44 

-0.1357 

0.0051 

0.0229 

171.81 

-4.34 

0.29 

-0.1357 

0.0051 

0.0254 

170.31 

-5.33 

9.07 

-0.1357 

0.0051 

0.0279 

168.85 

-6.83 

16.41 

-0.1357 

0.0051 

0.0305 

167.28 

-7.40 

23.07 

-0.1357 

0.0051 

0.033 

165.74 

-6.02 

28.96 

-0.1357 

0.0051 

0.0356 

164.37 

-5.97 

33.98 

-0.1357 

0.0051 

0.0381 

160.65 

-3.40 

42.66 

-0.1357 

0.0051 

0.0407 

157.26 

-26.34 

54.89 

-0.1357 

0.0051 

0.0432 

156.11 

-25.97 

55.07 

-0.1357 

0.0051 

0.0457 

156.86 

-24.85 

55.54 

-0.1357 

0.0051 

0.0483 

156.13 

-23.29 

55.59 

-0.1357 

0,0051 

0.0508 

156.31 

-21.36 

55.60 

-0.1357 

0.0025 

0.0076 

119.05 

52.27 

-33.48 

-0.1357 

0.0025 

0.0102 

176.95 

-33.16 

-43.47 

-0.1357 

0.0025 

0.0127 

175.83 

-32.26 

-35.45 

-0.1357 

0.0025 

0.0152 

174.82 

-7.26 

-25.71 

-0.1357 

0.0025 

0.0178 

173.73 

-2.53 

-16.27 

-0.1357 

0.0025 

0.0203 

172.45 

-0.03 

-7.07 

-0.1357 

0.0025 

0.0229 

171.31 

1.29 

1.45 

-0.1357 

0.0025 

0.0254 

169.91 

4.52 

9.03 

-0.1357 

0.0025 

0.0279 

168.41 

6.58 

16.06 

-0.1357 

0.0025 

0.0305 

166.68 

6.23 

22.21 

-0.1357 

0.0025 

0.033 

165.28 

7.75 

28.12 
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X 

y 

z 

u 

V 

w 

(m) 

(m) 

(m) 

(m/s) 

(m/s) 

(m/s) 

-0.1357 

0.0025 

0.0356 

164.07 

7.48 

32.97 

-0.1357 

0.0025 

0.0381 

161.27 

9.05 

40.39 

-0.1357 

0.0025 

0.0407 

157.44 

-34.57 

52.96 

-0.1357 

0.0025 

0.0432 

156.24 

-29.97 

53.89 

-0.1357 

0.0025 

0.0457 

155.53 

-22.86 

54.27 

-0.1357 

0.0025 

0.0483 

155.34 

-19.51 

54.50 

-0.1357 

0.0025 

0.0508 

153.78 

-21.40 

54.32 

-0.1357 

0 

0.0076 

138.88 

-21.78 

-20.43 

-0.1357 

0 

0.0102 

150.80 

-21.48 

-48.16 

-0.1357 

0 

0.0127 

150.39 

-21.20 

-39.13 

-0.1357 

0 

0.0152 

148.48 

-21.14 

-31.29 

-0.1357 

0 

0.0178 

173.03 

-36.22 

-12.91 

-0.1357 

0 

0.0203 

172.18 

-24.65 

-5.50 

-0.1357 

0 

0.0229 

170.68 

0.58 

2.72 

-0.1357 

0 

0.0254 

169.45 

4.87 

9.47 

-0.1357 

0 

0.0279 

168.04 

7.30 

15.92 

-0.1357 

0 

0.0305 

166.63 

-1.94 

22.13 

-0.1357 

0 

0.033 

165.27 

-1.12 

27.18 

-0.1357 

0 

0.0356 

163.13 

-4.51 

33.62 

-0.1357 

0 

0.0381 

159.36 

-16.06 

40.87 

-0.1357 

0 

0.0407 

158.14 

-41.54 

49.13 

-0.1357 

0 

0.0432 

155.03 

-32.55 

51.19 

-0.1357 

0 

0.0457 

154.11 

-23.27 

52.22 

-0.1357 

0 

0.0483 

153.46 

-20.37 

53.00 

-0.1357 

0 

0.0508 

153.36 

-24.14 

53.66 
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Table  A-2:  Experimental  Data  for  20"  Angle  of  Attack,  80%  Chord. 


X 

(m) 

y 

(m) 

z 

(m) 

u 

(m/s) 

V 

(m/s) 

w 

(m/s) 

0 

0.0711 

0.0051 

139.68 

18.47 

50.44 

0 

0.0711 

0.0102 

139.34 

18.28 

55.35 

0 

0.0711 

0.0152 

140.74 

16.04 

59.41 

0 

0.0711 

0.0203 

144.34 

9.47 

63.23 

0 

0.0711 

0.0254 

151.37 

-3.77 

67.42 

0 

0.0711 

0.0305 

151.90 

-7.54 

69.27 

0 

0.0711 

0.0356 

154.99 

-17.43 

75.01 

0 

0.0711 

0.0407 

158.40 

-25.42 

72.23 

0 

0.0711 

0.0457 

157.48 

-29.04 

77.22 

0 

0.0711 

0.0508 

154.98 

-21.65 

69.35 

0 

0.0711 

0.0559 

158.70 

-33.14 

73.72 

0 

0.0711 

0.061 

159.21 

-32.35 

71.23 

0 

0.0711 

0.0661 

155.45 

-13.40 

66.39 

0 

0.0711 

0.0762 

160.76 

-19.55 

64.40 

0 

0.0711 

0.0813 

160.26 

-18.36 

61.25 

0 

0.0711 

0.0864 

160.23 

-16.33 

60.08 

0 

0.0711 

0.0915 

159.92 

-15.54 

59.16 

0 

0.0711 

0.0965 

159.77 

-14.13 

58.40 

0 

0.0711 

0.1016 

160.59 

-14.59 

58.22 

0 

0.0661 

0.0152 

141.82 

13.69 

49.16 

0 

0.0661 

0.0203 

142.20 

10.06 

57.34 

0 

0.0661 

0.0254 

145.86 

4.99 

62.30 
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X 

y 

z 

u 

V 

w 

(m) 

(m) 

(m) 

(m/s) 

(m/s) 

(m/s) 

0 

0.0661 

0.0305 

155.87 

-20.17 

67.82 

0 

0.0661 

0.0356 

156.95 

-21.10 

70.44 

0 

0.0661 

0.0407 

161.83 

-33.05 

72.15 

0 

0.0661 

0.0457 

158.90 

-34.43 

79.55 

1  ^ 

0.0661 

0.0508 

156.44 

-25.83 

70.99 

1  ^ 

0.0661 

0.0559 

159.79 

-37.81 

74.02 

0 

0.0661 

0.061 

160.08 

-36.60 

71.14 

0 

0.0661 

0.0661 

156.08 

-16.29 

66.71 

0 

0.0661 

0.0711 

159.23 

-23.85 

65.76 

0 

0.0661 

0.0762 

161.28 

-22.11 

63.53 

0 

0.0661 

0.0813 

160.86 

-19.73 

60.29 

0 

0.0661 

0.0864 

162.02 

-21.19 

59.33 

0 

0.0661 

0.0915 

160.48 

-16.53 

57.85 

0 

0.0661 

0.0965 

161.85 

-18.51 

57.71 

0 

0.0661 

0.1016 

161.03 

-15.49 

57.17 

0 

0.061 

0.0254 

142.18 

9.23 

53.53 

0 

0.061 

0.0305 

145.64 

2.04 

60.87 

0 

0.061 

0.0356 

158.45 

-24.95 

69.02 

0 

0.061 

0.0407 

162.09 

-36.05 

70.08 

0 

0.061 

0.0457 

162.04 

-44.66 

79.55 

0 

0.061 

0.0508 

157.60 

-29.43 

71.94 

0 

0.061 

0.0559 

160.86 

-42.01 

74.09 

0 

0.061 

0.061 

161.17 

-39.97 

70.97  j 

0 

0.061 

0.0661 

156.52 

-18.52 

66.91  1 

0 

0.061 

0.0711 

159.75 

-25.89 

65.20  I 

0 

0.061 

0,0762 

161.85 

-23.79 

62.40  1 
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1  ^ 

y 

z 

u 

V 

w 

1 

(m) 

(m) 

(m/s) 

(m/s) 

(m/s) 

1  ^ 

0.061 

0.0813 

161.67 

-21.19 

58.75 

1  ^ 

0.061 

0.0864 

162.35 

-23.58 

57.84 

0 

0.061 

0.0915 

160.91 

-17.68 

56.57 

0 

0.061 

0.0965 

162.66 

-19.86 

56.65 

0 

0.061 

0.1016 

161.20 

-15.37 

55.86 

0 

0.0559 

0.0356 

146.32 

-3.02 

60.38 

0 

0.0559 

0.0407 

168.11 

-49.89 

70.08  1 

0 

0.0559 

0.0457 

164.30 

-50.83 

77.56 

0 

0.0559 

0.0508 

157.16 

-28.09 

68.71 

0 

0.0559 

0.0559 

163.18 

-47.49 

73.24 

0 

0.0559 

0.061 

162.35 

-44.63 

70.62 

0 

0.0559 

0.0661 

157.34 

-22.30 

67.12 

0 

0.0559 

0.0711 

161.03 

-29.80 

64.47 

0 

0.0559 

0.0762 

162.68 

-25.35 

60.93 

0 

0.0559 

0.0813 

162.79 

-22.82 

56.97 

0 

0.0559 

0.0864 

163.04 

-22.12 

56.09 

0 

0.0559 

0.0915 

161.45 

-17.89 

55.12 

0 

0.0559 

0.1016 

161.88 

-15.92 

54.68 

0 

0.0508 

0.0407 

171.23 

-59.15 

65.25  1 

0 

0.0508 

0.0457 

168.89 

-59.93 

68.61 

0 

0.0508 

0.0508 

161.30 

-38.48 

67.39 

0 

0.0508 

0.0559 

165.34 

-53.95 

70.95 

0 

0.0508 

0.061 

164.16 

-49.02 

68.53 

0 

0.0508 

0.0661 

158.53 

-26.47 

66.33 

0 

0.0508 

0.0711 

161.24 

-29.92 

62.99  I 

0 

0.0508 

0.0762 

163.12 

-25.61 

59.58  1 
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1  ^ 

(m) 

y 

(m) 

z 

(m) 

u 

(m/s) 

V 

(m/s) 

w  1 

(m/s)  1 

0 

0.0508 

0.0813 

163.65 

-23.39 

55.64  1 

0 

0.0508 

0.0864 

163.65 

-21.63 

54.02  1 

0 

0.0508 

0.0915 

161.88 

-17.10 

53.32 

0 

0.0508 

0.0965 

162.79 

-17.02 

53.62 

0.0508 

0.1016 

162.10 

-14.64 

53.45 

0 

0.0457 

0.0457 

168.89 

-69.85 

70.90 

0 

0.0457 

0.0508 

160.14 

-42.69 

64.41 

0 

0.0457 

0.0559 

166.56 

-58.79 

67.00 

0 

0.0457 

0.061 

165.39 

-51.62 

63.45 

0 

0.0457 

0.0661 

160.83 

-33.38 

64.33 

0 

0.0457 

0.0711 

163.05 

-32.08 

60.77 

0 

0.0457 

0.0762 

163.67 

-25.05 

57.14 

0 

0.0457 

0.0813 

164.61 

-22.37 

53.58 

0 

0.0457 

0.0864 

164.13 

-18.96 

51.84 

0 

0.0457 

0.0915 

162.52 

-14.71 

51.72 

0 

0.0457 

0.0965 

162.93 

-14.88 

52.05 

0 

0.0457 

0.1016 

162.69 

-13.12 

52.46 

0 

0.0407 

0.0457 

171.48 

-77.45 

66.14 

0 

0.0407 

0.0559 

166.38 

-60.13 

63.20 

0 

0.0407 

0.061 

165.74 

-52.61 

55.26 

0 

0.0407 

0.0661 

160.30 

-31.01 

58.18 

0 

0.0407 

0.0711 

163.55 

-33.70 

58.89 

0 

0.0407 

0.0762 

164.30 

-25.18 

54.87 

0 

0.0407 

0.0813 

165.46 

-22.54 

50.06 

0 

0.0407 

0.0864 

164.91 

-20.53 

49.97 

0 

0.0407 

0.0915 

163.28 

-14.74 

50.15 
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X 

y 

z 

u 

V 

w 

(m) 

(m) 

(m) 

(m/s) 

(tn/s) 

(m/s) 

0 

0.0407 

0.0965 

163.78 

-14.22 

50.69 

0 

0.0407 

0.1016 

162.10 

-12.57 

55.03 

0 

0.0356 

0.0559 

165.58 

-58.64 

50.92 

0 

0.0356 

0.061 

165.25 

-51.25 

48.56 

0 

0.0356 

0.0661 

161.97 

-32.18 

52.81 

0 

0.0356 

0.0711 

163.90 

-32.42 

54.99 

0 

0.0356 

0.0762 

164.32 

-23.28 

51.92 

0 

0.0356 

0.0813 

164.69 

-19.32 

47.58 

0 

0.0356 

0.0864 

165.33 

-21.45 

48.16 

0 

0.0356 

0.0915 

163.07 

-13.46 

48.30 

0 

0.0356 

0.0965 

163.61 

-13.60 

49.36 

0 

0.0356 

0.1016 

162.10 

-11.45 

53.15 

0 

0.0305 

0.0559 

166.61 

-55.23 

38.04 

0 

0.0305 

0.061 

165.61 

-47.63 

41.15 

0 

0.0305 

0.0661 

165.88 

-34.33 

46.80 

0 

0.0305 

0.0711 

164.85 

-29.14 

48.43 

0 

0.0305 

0.0762 

164.59 

-20.43 

48.88 

0 

0.0305 

0.0864 

165.77 

-18.04 

46.52 

0 

0.0305 

0.0915 

163.10 

-11.78 

47.01 

0 

0.0305 

0.0965 

163.24 

-11.35 

47.92 

0 

0.0305 

0.1016 

161.66 

-10.42 

52.76 

0 

0.0254 

0.0559 

166.85 

-48.65 

30.86 

0 

0.0254 

0.061 

166.25 

-42.74 

35.26 

0 

0.0254 

0.0661 

166.28 

-27.94 

41.48 

0 

0.0254 

0.0711 

165.43 

-23.59 

43.11 

0 

0.0254 

0.0762 

165.09 

-16.80 

_ 

44.97 
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X 

y 

z 

u 

V 

w 

(m) 

(m) 

(m) 

(m/s) 

(m/s) 

(m/s) 

0 

0.0254 

0.0813 

165.07 

-13.70 

42.91 

0 

0.0254 

0.0864 

165.11 

-13.55 

44.47 

0 

0.0254 

0.0915 

163.37 

-9.98 

45.54 

0 

0.0254 

0.0965 

155.27 

-9.44 

69.35 

0 

0.0254 

0.1016 

162.30 

-9.37 

50.52 

0 

0.0203 

0.0559 

167.35 

-42.03 

26.61 

0 

0.0203 

0.061 

166.52 

-36.87 

31.41 

0 

0.0203 

0.0661 

166.46 

-21.32 

36.99 

0 

0.0203 

0.0711 

165.51 

-16.68 

37.87 

0 

0.0203 

0.0762 

164.51 

-12.02 

43.80 

0 

0.0203 

0.0813 

165.39 

-13.45 

41.45 

0 

0.0203 

0.0864 

164.93 

-10.07 

43.33 

0 

0.0203 

0.0915 

163.59 

-8.78 

44.43 

0 

0.0203 

0,0965 

163.55 

-6.77 

45.54 

0 

0.0203 

0.1016 

161.99 

-6.73 

50.64 

0 

0.0152 

0.0457 

168.53 

-32.05 

11.85 

0 

0.0152 

0.0559 

167.78 

-34.98 

24.20 

0 

0.0152 

0.061 

167.01 

-31.00 

28.67 

0 

0.0152 

0.0661 

166.81 

-14.58 

33.34  1 

0 

0.0152 

0.0711 

165.55 

-17.59 

35.92 

0 

0.0152 

0.0762 

164.52 

-7.99 

40.91 

0 

0.0152 

0.0864 

164.43 

-7.65 

41.92 

0 

0.0152 

0.0915 

163.87 

-7.45 

43.33 

0 

0.0152 

0.0965 

163.58 

-7.27 

44.82 

1  ^ 

0.0152 

0.1016 

163.29 

-4.63 

46.32  1 

1  0 

0.0102 

0.0457 

168.39 

-29.46 

11.30  1 
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X 

y 

z 

u 

V 

w 

(m) 

(m) 

(m) 

(m/s) 

(m/s) 

(m/s) 

0 

0.0102 

0.0559 

167.75 

-27.29 

23.35 

0 

0.0102 

0.061 

167.05 

-24.00 

27.05 

0 

0.0102 

0.0661 

166.89 

-7.83 

31.37 

0 

0.0102 

0.0711 

165.87 

-6.60 

34.73 

0 

0.0102 

0.0762 

164.91 

-4.76 

37.99 

0 

0.0102 

0.0864 

164.46 

-0.79 

41.10 

0 

0.0102 

0.0915 

164.11 

-4.71 

42.68 

0 

0.0102 

0.0965 

163.48 

-5.82 

44.33 

0 

0.0102 

0.1016 

163.12 

-6.17 

45.86 

0 

0.0051 

0.0457 

168.80 

-21.05 

11.10 

0 

0.0051 

0.0559 

167.75 

-20.39 

22.20 

0 

0.0051 

0.061 

167.09 

-18.15 

26.54 

0 

0.0051 

0.0661 

167.07 

-1.18 

30.31  1 

0 

0.0051 

0.0711 

165.80 

-0.98 

34.74  1 

0 

0.0051 

0.0762 

164.62 

5.44 

36.50  1 

0 

0.0051 

0.0813 

164.52 

-0.37 

38.39 

0 

0.0051 

0.0864 

164.88 

-2.33 

40.64 

0 

0.0051 

0.0915 

164.18 

-2.04 

42.45 

0 

0.0051 

0.0965 

163.75 

-1.33 

44.22 

0 

0.0051 

0.1016 

162.86 

-2.36 

45.96 

0 

0 

0.0407 

168.66 

-12.39 

7.14 

0 

0 

0.0457 

168.45 

-13.48 

12.73 

0 

0 

0.0559 

167.86 

-13.27 

22.10 

0 

0 

0.061 

167.21 

-12.23 

26.15  1 

0 

0 

0.0661 

167.03 

5.01 

30.41  j 

0 

0 

0.0711 

165.95 

4.54 

34.53  1 
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X 

(m) 

y 

(m) 

z 

(m) 

u 

(m/s) 

V 

(m/s) 

w  1 

(m/s)  1 

0 

0 

0.0762 

165.15 

3.83 

35.84 

0 

0 

0.0813 

164.56 

2.22 

38.27 

0 

0 

0.0864 

164.73 

-0.99 

40.67  1 

1  ^ 

0 

0.0915 

163.96 

0.06 

42.72  1 

1  ^ 

0 

0.0965 

163.51 

-1.21 

44.32  1 

1  0 

0 

0.1016 

163.05 

-1.94 

45.57  1 
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Table  A-3:  Experimental  Data  for  30“  Angle  of  Attack,  40%  Chord. 


X 

y 

z 

u 

V 

w 

(m) 

(m) 

(m) 

(m/s) 

(m/s) 

(m/s) 

-0.1357 

0.0407 

0.0051 

116.04 

39.87 

63.53 

-0.1357 

0.0407 

0.0076 

117.59 

30.07 

68.27 

-0.1357 

0.0407 

0.0102 

118.22 

28.28 

72.39 

-0.1357 

0.0407 

0.0127 

113.2 

40.84 

74.99 

-0.1357 

0.0407 

0.0152 

115.41 

34.09 

79.99 

-0.1357 

0.0407 

0.0178 

116.81 

30.1 

83.86 

-0.1357 

0.0407 

0.0203 

119.86 

24.27 

86.86 

-0.1357 

0.0407 

119.71 

21.25 

92.41 

-0.1357 

0.0407 

0.0254 

121.38 

11.13 

92.38 

-0.1357 

0.0407 

0.0279 

122.58 

6.37 

93.71 

-0.1357 

0.0407 

0.0305 

128 

-4.41 

96.57 

-0.1357 

0.0407 

0.033 

129.2 

11.56 

96.49 

-0.1357 

0.0407 

0.0356 

130.05 

-14.11 

98.86 

-0.1357 

0.0407 

0.0381 

131.91 

-18.09 

99.07 

-0.1357 

0.0407 

0.0407 

132.69 

-20.61 

100 

-0.1357 

0.0407 

0.0432 

133.57 

-20.57 

99.74 

-0.1357 

0.0407 

0.0457 

134.46 

-23.58 

99.89 

-0.1357 

0.0407 

0.0483 

135.55 

-29.82 

100.51 

-0.1357 

0.0407 

0.0508 

135.9 

-27.15 

100.04 

-0.1357 

0.0381 

0.0051 

115.42 

20.27 

62.25 

-0.1357 

0.0381 

0.0076 

106.32 

48.2 

62.26 

-0.1357 

0.0381 

0.0102 

108.98 

40.48 

72.94 

-0.1357 

0.0381 

0.0127 

113.12 

27.31 

81.64 
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X 

y 

z 

u 

V 

w 

(m) 

(m) 

(m) 

(m/s) 

(m/s) 

(m/s) 

-0.1357 

0.0381 

0.0152 

117.3 

24.8 

85.34 

-0.1357 

0.0381 

0.0178 

117.93 

20.8 

89 

-0.1357 

0.0381 

0.0203 

120.56 

19.54 

90.19 

-0.1357 

0.0381 

0.0229 

120.32 

16.81 

94.06 

-0.1357 

0.0381 

0.0254 

122.05 

8.74 

94.73 

-0.1357 

0.0381 

0.0279 

126.7 

-6.65 

98.16 

1  -0.1357 

0.0381 

0.0305 

128.23 

-4.89 

99.31  1 

1  -0.1357 

0.0381 

0.033 

129.72 

-11.12 

99.72  ! 

-0.1357 

0.0381 

0.0356 

131.65 

-13.4 

100.2  1 

-0.1357 

0.0381 

0.0381 

131.95 

-7.44 

100.07  j 

-0.1357 

0.0381 

0.0407 

132.46 

-15.39 

100.25  1 

-0.1357 

0.0381 

0.0432 

133.24 

-17.91 

100.1 

-0.1357 

0.0381 

0.0457 

134.41 

-22.19 

100.12 

-0.1357 

0.0381 

0.0483 

134.9 

-26.98 

100.38  1 

-0.1357 

0.0381 

0.0508 

134.73 

-24.57 

98.48 

-0.1357 

0.0356 

0.0051 

99.37 

45.48 

46.53 

-0.1357 

0.0356 

0.0076 

108.62 

25.07 

62.82 

D 

-0.1357 

0.0356 

0.0102 

119.69 

15.71 

73.5 

-0.1357 

0.0356 

0.0127 

121.19 

10.9 

77.15 

-0.1357 

0.0356 

0.0152 

126.7 

-11.4 

82.69 

-0.1357 

0.0356 

0.0178 

124.31 

-1.26 

81.57 

-0.1357 

0.0356 

0.0203 

126.03 

-8.3 

87.48 

1  -0.1357 

0.0356 

0.0229 

127.28 

-5.2 

90.14 

-0.1357 

0.0356 

0.0254 

127,48 

10.18 

92.17 

-0.1357 

0.0356 

0.0279 

126.61 

4.54 

95.8 

-0.1357 

0.0356 

0.0305 

124.96 

-0.79 

98.69 
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X 

(m) 

y 

(m) 

u 

(m/s) 

V 

(m/s) 

w 

(m/s) 

-0.1357 

0.0356 

0.033 

126.31 

-5.33 

100.1 

-0.1357 

0.0356 

0.0356 

130.37 

-12.12 

101.54 

-0.1357 

0.0356 

0.0381 

131.45 

-17.48 

101.5 

-0.1357 

0.0356 

0.0407 

132.64 

-19.74 

101.52 

-0.1357 

0.0356 

0.0432 

133.77 

-22.83 

100.54 

-0.1357 

0.0356 

0.0457 

135.57 

-28.1 

100.54 

-0.1357 

0.0356 

0.0483 

136.17 

-29.06 

100.39 

-0.1357 

0.0356 

0.0508 

136.36 

-26.2 

99.5 

-0.1357 

0.033 

0.0076 

124.63 

19.5 

70.66 

-0.1357 

0.033 

0.0102 

124.25 

17.26 

73.5 

-0.1357 

0.033 

0.0127 

124.4 

20.46 

76.2 

-0.1357 

0.033 

0.0152 

124.49 

20.77 

78.8 

-0.1357 

0.033 

0.0178 

124.81 

15.35 

81.75 

-0.1357 

0.033 

0.0203 

121.05 

-16.49 

91.19 

-0.1357 

0.033 

0.0229 

125.09 

-22.74 

86.07 

-0.1357 

0.033 

0.0254 

125.72 

-19.8 

88.11 

-0.1357 

0.033 

0.0279 

124.34 

-15.22 

93.8 

-0.1357 

0.033 

0.0305 

125.26 

-7.2 

94.8 

-0.1357 

0.033 

0.033 

126.86 

-3.46 

97.29 

-0.1357 

0.033 

0.0356 

130.73 

-19.18 

102.03 

-0.1357 

0.033 

0.0381 

130.43 

-19.28 

102.76 

-0.1357 

0.033 

0.0407 

131.65 

-3.62 

102.9 

-0.1357 

0.033 

0.0432 

132.73 

-22.94 

101.44 

-0.1357 

0.033 

0.0457 

133.39 

-24.29 

101.91 

-0.1357 

0.033 

0.0483 

134.83 

-26.92 

101.31 

-0.1357 

0.033 

0.0508 

135.82 

-14.39 

100.65 
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X 

y 

z 

u 

V 

w 

(m) 

(m) 

(m) 

(m/s) 

(m/s) 

(m/s) 

-0.1357 

0.0305 

0.0102 

121.88 

-9.07 

70.42 

-0.1357 

0.0305 

0.0127 

122.86 

11.89 

73 

-0.1357 

0.0305 

0.0152 

123.88 

6.5 

75.75 

-0.1357 

0.0305 

0.0178 

125.11 

6.75 

78.81 

-0.1357 

0.0305 

0.0203 

125.76 

2.9 

81.27 

-0.1357 

0.0305 

0.0229 

126.26 

-1.02 

83.88 

-0.1357 

0.0305 

0.0254 

126.67 

-21.72 

85.79 

-0.1357 

0.0305 

0.0279 

126.99 

-31.32 

88.13 

-0.1357 

0.0305 

0.0305 

-25.13 

89.9 

-0.1357 

0.0305 

0.033 

126.93 

-17.4 

94.06 

-0.1357 

0.0305 

0.0356 

128.17 

-11.16 

95.91 

-0.1357 

0.0305 

0.0381 

129.8 

-16.16 

97.94 

-0.1357 

0.0305 

0.0407 

131.15 

-23.2 

100.7 

-0.1357 

0.0305 

0.0432 

130.91 

-23.56 

100.85 

-0.1357 

0.0305 

0.0457 

131.88 

-22.36 

100.77 

-0.1357 

0.0305 

0.0483 

132.83 

-24.66 

99.69 

-0.1357 

0.0305 

0.0508 

133.52 

-24.77 

98.7 

-0.1357 

0.0279 

0.0229 

125.32 

5.13 

80.25 

-0.1357 

0.0279 

0.0254 

126.26 

0.8 

82.7 

-0.1357 

0.0279 

0.0279 

127.41 

-14.44 

85.4 

-0.1357 

0.0279 

0.0305 

127.6 

-32.25 

87.6 

-0.1357 

0.0279 

0.033 

128.01 

-26.25 

90.14 

-0.1357 

0.0279 

0.0356 

132.8 

-29.11 

94.46 

-0.1357 

0.0279 

0.0381 

129.3 

-19.65 

^4.91 

-0.1357 

0.0279 

0.0407 

134.99 

-28.25 

99.58 

-0.1357 

0.0279 

0.0432 

135.41 

-34.31 

102.41 

66 


X 

y 

z 

u 

V 

w 

(m) 

(m) 

(m) 

(m/s) 

(m/s) 

(m/s) 

-0.1357 

0.0279 

0.0457 

135.06 

-33.83 

103.65 

-0.1357 

0.0279 

0.0483 

135.8 

-33.73 

103.45 

-0.1357 

0.0279 

0.0508 

136.67 

-18.87 

101.66  1 

-0.1357 

0.0254 

0.0279 

117.9 

3.36 

73.62 

-0.1357 

0.0254 

0.0305 

114.38 

-0.65 

76.97 

-0.1357 

0.0254 

0.033 

123.88 

-31.24 

84.41 

-0.1357 

0.0254 

0.0356 

126.93 

-30.24 

88.38 

1  -0.1357 

0.0254 

0.0381 

129.78 

-21.11 

92.22 

1  -0.1357 

0.0254 

0.0407 

139.09 

-40.87 

99.73  1 

1  -0.1357 

0.0254 

0.0432 

137.33 

-34.22 

99.8  1 

1  -0.1357 

0.0254 

0.0457 

137.47 

-36.37 

102.16  1 

1  -0.1357 

0.0254 

0.0483 

137.83 

-38.49 

102.72  1 

-0.1357 

0.0254 

0.0508 

137.58 

-38.73 

101.98 

-0.1357 

0.0229 

0.0381 

146.85 

-67.73 

109.62 

-0.1357 

0.0229 

0.0407 

150.96 

-58.47 

100.16 

-0.1357 

0.0229 

0.0432 

136.13 

-41.37 

96.49 

-0.1357 

0.0229 

0.0457 

136.23 

-37.55 

98.04 

-0.1357 

0.0229 

0.0483 

137.15 

-43.28 

99.88 

-0.1357 

0.0229 

0.0508 

138.43 

-46.08 

99.61 

-0.1357 

0.0203 

0.0356 

125.08 

-33.24 

91.83 

-0.1357 

0.0203 

0.0381 

146.39 

-77.52 

111.96 

-0.1357 

0.0203 

0.0407 

148.4 

-73.98 

107.04 

1  -0.1357 

0.0203 

0.0432 

149.99 

-59.3 

100.26 

1  -0.1357 

0.0203 

0.0457 

146.1 

-51.65 

104.15 

1  -0.1357 

0.0203 

0.0483 

144.32 

-44.29 

100.44 

1  -0.1357 

0.0203 

0.0508 

137.76 

-24.68 

101.3 
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X 

y 

z 

u 

V 

w 

(m) 

(m) 

(m) 

(m/s) 

(m/s) 

(m/s) 

-0.1357 

0.0178 

0.0381 

128.66 

-46.27 

95.66 

-0,1357 

0.0178 

0.0407 

148.73 

-80.15 

106.63 

-0.1357 

0.0178 

0.0432 

149.66 

-73.27 

101.46 

-0.1357 

0.0178 

0.0457 

147.2 

-53.9 

95.04 

-0.1357 

0.0178 

0.0483 

146.27 

-41.22 

92.2 

-0.1357 

0.0178 

0.0508 

147 

-39. 14 

87.26 

-0.1357 

0.0152 

0.0407 

145.45 

-64.45 

87.07 

-0.1357 

0.0152 

0.0432 

154.21 

-65.02 

81.02 

-0.1357 

0.0152 

0.0457 

155.01 

-56.19 

73.87 

-0.1357 

0.0152 

0.0483 

156.47 

-45.4 

67.39 

-0.1357 

0.0152 

0.0508 

153.5 

-37.85 

72.45 

-0.1357 

0.0127 

0.0407 

140.72 

-32.91 

52.09 

-0.1357 

0.0127 

0.0432 

158.32 

-55.65 

61.34 

-0.1357 

0.0127 

0.0457 

157.23 

-49.3 

62.95 

-0.1357 

0,0127 

0.0483 

157.43 

-45.54 

60.31 

-0.1357 

0.0127 

0.0508 

155.87 

-35.82 

63.83 

-0.1357 

0.0102 

0.0407 

140.24 

-26.26 

37.77 

-0.1357 

0.0102 

0.0432 

158.64 

-50.68 

52.3 

-0.1357 

0.0102 

0.0457 

158.37 

-43.31 

54.41 

-0.1357 

0.0102 

0.0483 

158.3 

-49.83 

55.88 

-0.1357 

0.0102 

0.0508 

156.89 

-37.41 

58.41 

-0.1357 

0.0076 

0.0432 

160.7 

-57.15 

43.02 

-0.1357 

0.0076 

0.0457 

160.06 

-44.33 

47.43 

-0.1357 

0.0076 

0.0483 

159.12 

-44.15 

50.92 

-0.1357 

0.0076 

0.0508 

157.65 

-36.04 

54.4 

-0.1357 

0.0051 

0.0432 

163.04 

-49.54 

36.03 
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X 

y 

z 

u 

V 

w 

(m) 

(m) 

(m) 

(m/s) 

(m/s) 

(m/s) 

-0.1357 

0.0051 

0.0457 

161.87 

-44.5 

41.46 

-0.1357 

0.0051 

0.0483 

159.93 

-45.04 

46.84 

-0.1357 

0.0051 

0.0508 

158.39 

-39.03 

51.18 

-0.1357 

0.0025 

0.0432 

164.69 

-48.63 

30.68 

-0.1357 

0.0025 

0.0457 

162.55 

-41.92 

38.21 

-0.1357 

0.0025 

0.0483 

160.65 

-40.22 

44.37 

-0.1357 

0.0025 

0.0508 

159.07 

-33.67 

49.16 

-0.1357 

0 

0.0407 

166.98 

-48.23 

21.28 

-0.1357 

0 

0.0432 

164.67 

-41.81 

29.72 

-0.1357 

0 

0.0457 

162.91 

-38.05 

37.27 

-0.1357 

0 

0.0483 

161.12 

-37.29 

43.24 

-0.1357 

0 

0.0508 

159.7 

-26.76 

47.87 

Table  A-4:  Experimental  Data  for  30“  Angle  of  Attack,  80%  Chord. 


X 

(m) 

y 

(m) 

z 

(m) 

u 

(m/s) 

V 

(m/s) 

w 

(m/s) 

0 

0.0711 

0.0051 

80.79 

92.56 

36.09 

0 

0.0711 

0.0102 

67.07 

132.87 

36.15 

0 

0.0711 

0.0152 

87.22 

92.14 

53.44 

0 

0.0711 

0.0203 

113.63 

29.9 

72.35 

0 

0.0711 

0.0254 

110.57 

24.03 

85.26 

0 

0.0711 

0.0305 

143.24 

-27.83 

114.45 

0 

0.0711 

0.0356 

128.33 

-3.24 

116.05 

0 

0.0711 

0.0407 

122.43 

32.34 

90.19 

0 

0.0711 

0.0457 

124.04 

16.76 

93.64 

0 

0.0711 

0.0508 

129.1 

2.66 

98.33 

0 

0.0711 

0.0559 

130.4 

-4.17 

100.11 

0 

0.0711 

0.061 

130.55 

-3.89 

102.02 

0 

0.0711 

0.0661 

132.6 

-10.48 

103.57 

0 

0.0711 

0.0711 

133.57 

-13.2 

102.77 

0 

0.0711 

0.0762 

133.94 

-19.26 

102.45 

0 

0.0711 

0.0813 

134.14 

-13.73 

101.02 

0 

0.0711 

0.0864 

136.27 

-20.62 

100.14 

0 

0.071 1 

0.0915 

137.24 

-19.27 

100.03 

0 

0.0711 

0.0965 

138.51 

-25.87 

99.55 

0 

0.0711 

0.1016 

139.51 

-27.54 

98.78 

0 

0.0661 

0.0305 

121.09 

5.29 

66.99 

0 

0.0661 

0.0356 

120.45 

3.45 

74.67 

0 

0.0661 

0.0407 

121.48 

-3.57 

82.8 
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X 

(m) 

y 

(m) 

z 

(m) 

u 

(m/s) 

V 

(m/s) 

w 

(m/s) 

0 

0.0661 

0.0457 

124.55 

-6.72 

87.9 

0 

0.0661 

0.0508 

126.31 

-2.35 

91.59 

0 

0.0661 

0.0559 

137.52 

-38.68 

98.9 

0 

0.0661 

0.0661 

135.41 

-37.76 

100.18 

0 

0.0661 

0.0711 

136.34 

-35.34 

100.56 

0 

0.0661 

0.0762 

137.09 

-37.07 

100.76 

0 

0.0661 

0.0813 

141.96 

-43.5 

93.91 

0 

0.0661 

0.0864 

139.09 

-48.  i4 

101.7 

0 

0.0661 

0.0915 

139.65 

-47.05 

102.22 

0 

0.0661 

0.0965 

142.42 

-53.02 

101.94 

0 

0.0661 

0.1016 

130.11 

-24.69 

92.84 

0 

0.061 

0.0457 

123.45 

16.63 

78.97 

0 

0.061 

0.0508 

125.08 

13.95 

85.59 

0 

0.061 

0.0559 

127.56 

7.67 

89.69 

0 

0.061 

0.061 

128.96 

2.84 

92.58 

0 

0.061 

0.0661 

129.3 

-1.45 

94.71 

0 

0.061 

0.0711 

130.05 

-4.23 

96.2 

0 

0.061 

0.0762 

130.96 

-7.13 

96.96 

0 

0.061 

0.0813 

137.13 

-39.71 

100.33 

0 

0.061 

0.0864 

137.14 

-45.09 

100.27 

0 

0.061 

0.0915 

137.66 

-43.53 

100.48 

0 

0.061 

0.0965 

138.86 

-45.76 

99.61 

0 

0.061 

0.1016 

141.1 

-46.17 

99.17 

0 

0.0559 

0.0559 

132.19 

-11.44 

84.8 

0 

0.0559 

0.061 

134.2 

-14.08 

90.45 

0 

0.0559 

0.0661 

137.37 

-19.49 

95.41 
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X 

y 

z 

u 

V 

w 

(m) 

(m) 

(m) 

(tn/s) 

(m/s) 

(m/s) 

0 

0.0559 

0.0711 

142.64 

-36.71 

100.9 

0 

0.0559 

0.0762 

141.85 

-33.89 

101.66 

0 

0.0559 

0.0813 

136.41 

-22.43 

99.37 

0 

0.0559 

0.0864 

138.64 

-26.97 

100.89 

0 

0.0559 

0.0915 

138.51 

-29.23 

100.78 

0 

0.0559 

0.0965 

139.63 

-30.4 

100.12 

0 

0.0559 

0.1016 

140.64 

-31.4 

98.58 

0 

0.0508 

0.061 

157.12 

-74.34 

94.46 

0 

0.0508 

0.0661 

126.95 

-0.33 

83.38 

0 

0.0508 

0.0711 

128.97 

-4.15 

88.4 

0 

0.0508 

0.0762 

131.45 

-14.88 

92.57 

0 

0.0508 

0.0813 

132.37 

-18.19 

95.24 

0 

0.0508 

0.0864 

133.96 

-17.55 

95.13 

0 

0.0508 

0.0915 

134.4 

-22.07 

97.28 

0 

0.0508 

0.0965 

133.66 

-18.1 

95.72 

0 

0.0508 

0.1016 

135.43 

-21.33 

93.71 

0 

0.0457 

0.0661 

145.14 

-73.89 

114.39 

0 

0.0457 

0.0711 

155.66 

-62.16 

96.03 

0 

0.0457 

0.0762 

129.82 

-8.53 

86.54 

0 

0.0457 

0.0813 

131.47 

-14.83 

91.63 

0 

0.0457 

0.0864 

136.92 

-26.57 

94.63 

0 

0.0457 

0.0915 

136.48 

-27.56 

97.25 

0 

0.0457 

0.1016 

141.13 

-34.12 

95.21 

0 

0.0407 

0.0661 

146.3 

-93.39 

113.16 

0 

0.0407 

0.0711 

145.63 

-84.42 

112.69 

0 

0.0407 

0.0762 

145.29 

-71.’ 

112.48 
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X 

(m) 

y 

(m) 

z 

(m) 

u 

(m/s) 

V 

(m/s) 

[=“== 

w 

(m/s) 

0 

0.0407 

0.0813 

145.47 

-63.82 

111.55 

0 

0.0407 

0.0864 

138.69 

-34.45 

93.53 

0 

0.0407 

0.0915 

136.99 

-31.77 

94.69 

0 

0.0407 

0.0965 

139.76 

-34.67 

95.07 

0 

0.0407 

0.1016 

142.61 

-47.02 

94.75 

0 

0.0356 

0.0711 

146.49 

-92.58 

111.32 

0 

0.0356 

0.0762 

145.96 

-86.61 

1 10.92 

0 

0.0356 

0.0813 

140.22 

-57.04 

107.27 

0 

0.0356 

0.0864 

145.46 

-47.92 

92.35 

0 

0.0356 

0.0915 

139.63 

-33.32 

90.49 

0 

0.0356 

0.0965 

136.89 

-27.19 

90.55 

0 

0.0356 

0.1016 

139.41 

-29.65 

90.95 

0 

0.0305 

0.0762 

141.66 

-66.29 

92.16 

0 

0.0305 

0.0813 

147.62 

-40.54 

76.77 

0 

0.0305 

0.0864 

146.91 

-38.17 

78.47 

0 

0.0305 

0.0915 

147.17 

-36.91 

78.87 

0 

0.0305 

0.0965 

146.13 

-37.64 

82.23 

0 

0.0305 

0.1016 

143.67 

-32.58 

84.54 

0 

0.0254 

0.0762 

149.83 

-53.54 

63.46 

0 

0.0254 

0.0813 

150.9 

-41.9 

62.93 

0 

0.0254 

0.0864 

147.5 

-28.75 

66.11 

0 

0.0254 

0.0915 

150.1 

-32.41 

67.45  1 

0 

0.0254 

0.0965 

149.86 

-34.68 

70.8  I 

0 

0.0254 

0.1016 

149.78 

-30.78 

71.27  1 

0 

0.0203 

0.0813 

148.66 

-33.18 

55.3  1 

0 

0.0203 

0.0864 

152.33 

-35.68 

55.84  1 
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1  X 

1 

y 

(m) 

z 

(m) 

0 

0.0203 

0.0915 

0 

0.0203 

0.0965 

0 

0.0203 

0.1016 

1  ^ 

0.0152 

0.0762 

1  ^ 

0.0152 

0.0813 

1  ^ 

0.0152 

0.0864 

1  ^ 

0.0152 

0.0915 

1  ^ 

0.0152 

0.0965 

0 

0.0152 

0.1016 

0 

0.0102 

0.0711 

0 

0.0102 

0.0762 

0 

0.0102 

0.0813 

0 

0.0102 

0.0864 

0 

0.0102 

0.0915 

0 

0.0102 

0.0965 

0 

0.0102 

0.1016 

0.0051 

0.0661 

0 

0.0051 

0.0711 

0 

0.0051 

0.0762 

0 

0.0051 

0.0813 

0 

0.0051 

0.0864 

0 

0.0051 

0.0915 

1  ^ 

0.0051 

0.0965 

1  0 

0.0051 

0.1016 

1  ^ 

0 

0.0457 

I  0 

0 

0.0508 

u 

(m/s) 


154.95 


156.98 


158.75 


139.45 


143.71 


148.28 


151.48 


153.67 


155.39 


123.48 


131.84 


137.16 


151.3 


153.67 


111.22 


119.59 


127.27 


133.1 


139.56 


144.47 


149.02 


152.34 


81.98 


88.3 


V 

(m/s) 


-32.75 


-33.57 


-30.16 


-51.99 


-43.52 


-39.4 


-30.19 


-31.34 


-28.66 


-56.38 


-45.4 


-42.56 


8.85 


.88 


-30.66 


-31.98 


-49.41 


-40.62 


-41.48 


-43. 14 


-37.9 


-29.95 


-21.52 


-25.76 


-35.16 


-36.76 


w 

(m/s) 


55.54 


54.33 


52.72 


68.54 


64.64 


60.82 


59.39 


57.55 


56.34 


87.3 


79.67 


74.35 


69.42 


65.01 


60.54 


58.43 


102 


93.85 


86.37 


80.95 


74.42 


68.7 


63.36 


60.27 


130.4 


124.87 
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X 

(m) 

y 

(m) 

z 

(m) 

u 

(m/s) 

V 

(m/s) 

w 

(m/s) 

0 

0 

0.0559 

95.11 

-35.87 

118.73 

0 

0 

0.061 

102.13 

-34.74 

111.45  1 

0 

0 

0.0661 

109.88 

-33.25 

104.38 

0 

0 

0.0711 

117.59 

-36.07 

97.21 

0 

0 

0.0762 

124.94 

-40.09 

90.31 

0 

0 

0.0813 

132.08 

-31.22 

82.49 

0 

0 

0.0864 

138.49 

-23.71 

76.48 

0 

0 

0.0915 

143.62 

-19.28 

71.26 

0 

0 

0.0965 

148.52 

-16.18 

66.04 

0 

0 

0.1016 

152.08 

-14.71 

62.3 
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Table  A-5:  Experimental  Data  for  35“  Angle  of  Attack,  40%  Chord. 


X 

y 

z 

u 

V 

w 

(m) 

(m) 

(m) 

(m/s) 

(m/s) 

(m/s) 

-0.1357 

0.0661 

0.0051 

112.15 

5.15 

77.06 

-0.1357 

0.0661 

0.0076 

112.6 

4.22 

78.2 

-0.1357 

0.0661 

0.0102 

112.74 

3.16 

79.18 

-0.1357 

0.0661 

0.0127 

113.38 

4.54 

80.41 

-0.1357 

0.0661 

0.0152 

114.38 

3.08 

81.99 

-0.1357 

0.0661 

0.0178 

115.49 

-3.91 

83.64 

-0.1357 

0.0661 

0.0203 

115.99 

-3.97 

84.79 

-0.1357 

0.0661 

0.0229 

116.75 

-8.87 

86.2 

-0.1357 

0.0661 

0.0254 

117.38 

-10.09 

87.5 

-0.1357 

0.0661 

0.0279 

118.19 

-13.93 

88.57 

-0.1357 

0.0661 

0.0305 

118.19 

-9.92 

89.55 

-0.1357 

0.0661 

0.033 

119.22 

-12.8 

90.58 

-0.1357 

0.0661 

0.0356 

119.83 

-9.8 

91.87 

-0.1357 

0.0661 

0.0381 

120.1 

-16.29 

92.56 

-0.1357 

0.0661 

0.0407 

120.36 

-19.93 

93.36 

-0.1357 

0.0661 

0.0432 

120.75 

-19.39 

94.53 

-0.1357 

0.0661 

0.0457 

121.14 

-20.24 

95.51 

-0.1357 

0.0661 

0.0483 

121.6 

-20.44 

96.42 

-0.1357 

0.0661 

0.0508 

122.05 

-21.38 

97.42 

-0.1357 

0.0661 

0.0534 

94.51 

30.71 

85.24 

-0.1357 
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-0.1357 

0.0279 

0.0305 

100.45 

13.13 

76.5 

-0.1357 

0.0279 

0.033 

107.01 

-24.56 

104.94 

-0.1357 

0.0279 

0.0356 

107.95 

-24.31 

108.77 

-0.1357 

0.0279 

0.0381 

109.77 

-14.97 

103.78 

-0.1357 

0.0279 

0.0407 

121.7 

-32.78 

107.57 

-0.1357 

0.0279 

0.0432 

121.66 

-43.13 

119.66 

-0.1357 

0.0279 

0.0483 

124.97 

-46.37 

111.41 

-0.1357 

0.0279 

0.0508 

121.54 

-42.86 

110.47 

-0.1357 

0.0279 

0.0534 

98.46 

9.19 

95.51 

-0.1357 

0.0279 

0.0559 

102.04 

3.35 

95.88  1 

-0.1357 

0.0279 

0.0584 

97.82 

17.19 

86.38  1 

-0.1357 

0.0279 

0.061 

107.73 

7.8 

83.73 

-0.1357 

0.0279 

0.0635 

112.89 

0.52 

81.94 

-0.1357 

0.0254 

0.0356 

94.37 

5.14 

92.11 

-0.1357 

0.0254 

0.0381 

97.91 

-10.13 

103.96 

-0.1357 

0.0254 

0.0407 

97.43 

10.88 

92.69  1 

-0.1357 

0.0254 

0.0432 

100.76 

10.77 

95.28  1 

-0.1357 

0.0254 

0.0457 

111.28 

-16.73 

107.84  1 

-0.1357 

0.0254 

0.0483 

108.66 

-11.47 

103.75  1 

-0.1357 

0.0254 

0.0508 

106.54 

-12.04 

100.46  1 
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X 

y 

z 

u 

V 

w 

(m) 

(m) 

(m) 

(m/s) 

(m/s) 

(m/s) 

-0.1357 

0.0254 

0.0534 

96.07 

12.57 

91.69 

-0.1357 

0.0254 

0.0559  ! 

106.91 

1.65 

90.09 

-0.1357 

0.0254 

0.0584 

100.92 

18.46 

79.99 

-0.1357 

0.0254 

0.061 

109.85 

8.32 

79.67 

-0.1357 

0.0254 

0.0635 

116.88 

-3.65 

80.16 

-0.1357 

0.0229 

0.0432 

94.58 

14.8 

88.4 

-0.1357 

0.0229 

0.0457 

110.13 

-4,17 

89.08 

-9.1357 

0.0229 

0.0483 

114.86 

-17.38 

101.41 

-0.1357 

0.0229 

0.0508 

116.29 

-32.28 

106.74 

-0.1357 

0.0229 

0.0534 

111.98 

-6.19 

91.81 

-0.1357 

0.0229 

0.0559 

114.82 

-9.65 

85.61 

-0.1357 

0.0229 

0.0584 

108.12 

9.28 

78.39 

-0.1357 

0.0229 

0.061 

120.81 

-P  23 

82.51 

-0.1357 

0.0229 

0.0635 

124.54 

-13.13 

80.69 

-0.1357 

0.0203 

0.0432 

79.31 

29.09 

83.98 

-0.1357 

0.0203 

0.0457 

77.63 

34.65 

82.08 

-0.1357 

0.0203 

0.0483 

81.16 

47.47 

73.72 

-0.1357 

0.0203 

0.0508 

114.69 

-15.83 

96.19 

-0.1357 

0.0203 

0.0534 

95.43 

25.92 

80.76 

-0.1357 

0.0203 

0.0559 

96.48 

23.46 

75.4 

-0.1357 

0.0203 

0.0584 

109.58 

12.47 

71.99 

-0.1357 

0.0203 

0.061 

122.32 

-11.39 

77.69 

-0.1357 

0.0203 

0.0635 

123.07 

-15.42 

78.16 

-0.1357 

0.0178 

0.0457 

112.54 

-24.66 

85.63 

-0.1.57 

0.0178 

0.0483 

i!3  53 

-12.8 

76.54 

-0.1357 

0.0178 

0.0508 

115.92 

L 

1.27 

70.17 
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X 

(m) 

y 

(m) 

z 

(m) 

u 

(m/s) 

V 

(m/s) 

w 

(m/s) 

-0.1357 

0.0178 

0.0534 

115.2 

8.22 

69.06 

-0.1357 

0.0178 

0.0559 

118.11 

-2.61 

73.13 

-0.1357 

0.0178 

0.0584 

104.83 

27.58 

63.21 

-0.1357 

0.0178 

0.061 

118.48 

0.81 

69.78 

-0.1357 

0.0178 

0.0635 

123.71 

-7.69 

73.09 

-0.1357 

0.0152 

0.0457 

104.09 

16.77 

56.82 

-0.1357 

0.0152 

0.0483 

108.46 

12.14 

57.28 

-0.1357 

0.0152 

0.0508 

111.07 

12.55 

56.61 

-0.1357 

0.0152 

0.0534 

102.81 

33.27 

51.32 

-0.1357 

0.0152 

0.0559 

105.98 

28.8 

56.92 

-0.1357 

0.0152 

0.0584 

104.23 

31.11 

57.08 

-0.1357 

0.0152 

0.061 

109.28 

20.3 

59.88 

-0.1357 

0.0152 

0.0635 

117.97 

2.83 

65.38 

-0.1357 

0.0127 

0.0483 

109.6 

22.4 

45.49 

-0.1357 

0.0127 

0.0508 

117.9 

4.16 

52.29 

-0.1357 

0.0127 

0.0534 

105.71 

28.32 

46.23 

-0.1357 

0.0127 

0.0584 

106.32 

35.41 

50.89 

-0.1357 

0.0127 

0.061 

105.31 

33.96 

51.71 

-0.1357 

0.0127 

0.0635 

112.49 

19.4 

57.05 

-0.1357 

0.0102 

0.0483 

117.99 

10.02 

41.16 

-0.1357 

0.0102 

0.0508 

120.37 

9.96 

45.49 

-0.1357 

0.0102 

0.0534 

117.29 

11.18 

46.91 

-0.1357 

0.0102 

0.0559 

113.45 

20.43 

46.13 

-0.1357 

0.0102 

0.0584 

112.64 

26.02 

47.87 

-0.1357 

0.0102 

0.061 

113 

25.97 

50.61 

-0.1357 

0.0102 

0.0635 

113.6 

22.37 

53.55 
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X 

y 

z 

u 

V 

w 

(m) 

(m) 

(m) 

(m/s) 

(m/s) 

(m/s) 

-0.1357 

0.0076 

0.0508 

122.15 

12.13 

35.99 

-0.1357 

0.0076 

0.0559 

114.87 

21.27 

40.17 

-0.1357 

0.0076 

0.0584 

115.76 

20.92 

43.81 

-0.1357 

0.0076 

0.061 

112.66 

26.22 

45.26 

-0.1357 

0.0076 

0.0635 

111.04 

33.54 

47.35 

-0.1357 

0.0051 

0.0483 

129.59 

4.96 

24.07 

-0.1357 

0.0051 

0.0508 

130.61 

7.1 

31.07 

-0.1357 

0.0051 

0.0534 

127.9 

11.11 

35.63 

-0.1357 

0.0051 

0.0559 

124.41 

15.44 

38.85 

-0.1357 

0.0051 

0.0584 

120.03 

18.36 

41.01 

-0.1357 

0.0051 

0.061 

116.95 

26.35 

43.04 

-0.1357 

0.0051 

0.0635 

123.14 

13.82 

51.36 

-0.1357 

0.0025 

0.0508 

125.93 

24.42 

20.43 

-0.1357 

0.0025 

0.0534 

124.55 

28 

26.24 

-0.1357 

0.0025 

0.0559 

120.32 

31.07 

30.61 

-0.1357 

0.0025 

0.0584 

119,73 

26.34 

36.77 

-0.1357 

0.0025 

0.061 

115.02 

33.55 

39.04 

-0.1357 

0.0025 

0.0635 

115.38 

31.06 

44.07 

-0.1357 

0 

0.0635 

104.5 

71.6 

36.93 
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Laser  Velocimetry  Seed  Particles  Within 
Compressible,  Vortical  Flows 

Mark  S.  Maurice* 

Wright  Laboratory,  Wright-Paiterson  Air  Force  Base,  Ohio  45433 


The  ablUly  of  (ccd  partkla  to  penetrate  and  accarateiy  track  vonkce  la  of  critical  Importance  to  Ihc 
analyrit  of  laecr  tehnlmttry  (LV)  mcasnrenicnti  wMtan  tbeae  Bow  rimetnret.  la  applyint  a  perllcte 
eqnatloo  of  motion  to  vortM  flowfleMi  which  extend  Into  the  mpereonlc,  compreerible  rethne.  two 
approachci  arc  contldcrcd.  fint.  an  Ideal,  potentlai  vortex  lx  develop  for  a  comprearihle  flowfield.  Ax 
an  aid  to  the  dcxlgn  and  aaalyxix  of  vottkal  llowflcid  xnrvcyx  over  a  wide  ranpe  of  Independent  param 
cterx,  Chix  model  lx  axed  to  plol  the  time  aad  poxitlott  at  which  any  partldc  xtartint  at  reit  wMhla  a  vortex 
win  bcgla  to  track  the  flow  velocity  wHhla  3%  error.  As  a  xpM^  appUcatloa,  the  potential  vortex  lx 
then  axed  to  estimate  the  dynamic  btoa  of  LV  mcasaremcntx  taken  wHhla  a  vortex  riied  hnm  a  7S>dcg 
dcMa  wli«  at  M  den  •■tic  of  attack  In  a  Mach  1.9  supersonic  flow.  The  second  approach  asm  a  com- 
pntoliooally  derived  Navter-Stokcx  flowfleld  oohilloo  la  place  at  the  potential  vortex  atodcL  The  com- 
pnlalfonal  flowfleld  method  pcadletx  that  acevrate  LV  mcamrcmfatx  wHhte  the  delta  whin  flowfldd 
qaire  seed  partlekx  no  larncr  than  t.l-t.Z  ftm  la  diameter,  and  deflnex  the  anxeeded  Inner  vortex  caca 
rcplen.  Bath  approaches  show  an  Inertaxc  la  vetocUy  bias  which  lx  nearly  proportloani  to  particle  dl- 
mnetcr,  xtrcxxlan  the  need  for  a  monodlxpcrac  seed  of  known  xixe  to  resolve  parllde  btos  In  rnmplix 
flows. 


Nomcnclatan 

Cd  w  drag  caefTicient.  Eq.  (IS) 

C<f|  =  function  defined  by  Eq.  (16) 

Cp  >=  specific  heat 

d  °  putkle  diameter 

C  =  function  defined  by  Eq.  (17) 

h  =■  function  defined  by  Eq.  (18) 

L  ”  root  chord  length 

M  »  Mach  number 

m  •  matt 

P  •  pressure 

R  =  gas  constant 

/fro  ~  Reynolds  number,  Eq.  (2) 

r  -  radial  special  coordinate 

f  =  transformed,  dimensionless  radial  coordinate.  Eq. 

(12) 

5  >  Sutherland  viscosity  constant 

St  >  Stokes  number.  Eq.  (4) 

t  >  distance,  used  in  (6) 

T  =  temperature 

r  “  time 

u  fluid  velocity 

V  =  particle  velocity 

w  w  modified  panicle  relative  velocity,  Eq.  (14) 

x,y,t  •  Cartesian  spacial  coordinates 

%ERR  -  measure  of  panicle  velocity  bias,  Eqs.  (24)  and 
(25) 

a  -  angle  of  attack 

B  »  density  ratio,  Eq.  (3) 

r  ~  circulation 

y  ratio  of  specific  heats 

0  “  angular  spacial  coordinate  in  a  cylindrical  system 

p  •  molecular  viscosity 

ir  »  dimensionless  variable  groups.  Eqs.  (26-28) 
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p  •  density 

r  =  dimensionless  time.  Eq.  (23) 

ai  vottkity 

Subteriptt 

f  •  fluid  property 

ij  w  Cartesian  coordinate  indkial  notation 

o  w  propeny  in  the  outer  field  of  a  two-dimensional 
vortex 

p  =  paitkie  property 

R  »  the  relative  difference  between  the  pnttkie  and  the 

fluid 

0  >  initial  condition 


InlroductitNi 

HE  dynamks  of  small  pankles  within  aerodynamk  flow- 
fklds  are  of  interest  for  a  variety  of  applkationt.  In  pro¬ 
pulsive  engines,  for  example,  combustion  partkies  can  have 
an  erosive  effect  on  turbine  blades,  as  well  as  contributing  to 
performance  losses  through  the  exit  noezle.  External  aerody- 
namk  and  control  surface  degradation  can  be  strong  fiinctiont 
of  almospherk  dispersion,  and  sediment  Irampon  can  be  mod¬ 
eled  by  the  dynamks  of  individual  partkies.  In  the  application 
of  flow  visu^izalion  or  laser  vehximeliy,  however,  it  is  re¬ 
quired  that  the  pattkle  dynamics  of  die  s^  material  rlescribe 
the  fluid  dynamks  of  the  flowfleld.  The  seed  partkies  are  as¬ 
sumed  to  follow  the  fluid  streamlines,  and  not  to  alter  the  cor¬ 
responding  unseeded  flow. 

In  flow  legiont  with  large  gradienis  such  at  near  shock  waves, 
or  regioru  with  high  turbulence  frequencies,  pnrtkie  dynamk 
bias  can  dominate  laser  vehximeliy  (LV)  measurements.  For 
the  case  of  a  seed  paitkie  passing  through  a  shock  wave,  it 
has  been  shown  by  Nkhols'  that  the  paitkie  lag  can  be  sub¬ 
stantial.  if  the  mean  velocity  and  turbulence  intensity  are  de¬ 
termined  from  ensembled  averages  from  a  variety  of  particle 
sizes  passing  through  the  shock  wave  and  measured  at  a  single 
downstream  location,  these  quantities  may  be  grossly  over¬ 
estimated.  On  the  other  hand,  it  has  been  shown  by  Hjelmfelt 
and  Mockios'  that  a  particle  subjected  to  a  high-frequency, 
oscillating  flowfleld  may  greatly  underestimate  the  tuiteleixc 
intensity  of  the  fluid. 
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APPENDIX  B:  PARTICLE  DYNAMICS  IN  VORTICAL  FLOWS 


This  appendix  contains  details  of  the  particle  dynamics  in  vortical  flows.  It  is  a 
reproduction  of  an  AIAA  Journal  article  (Ref  11)  written  as  a  work  of  the  Aeromechanics 
Division,  Flight  Dynamics  Directorate,  Wright  Laboratory,  Wright-Patterson  Air  Force  Base. 
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Allhough  fluid  streamlines  are  cireular  within  a  two-di¬ 
mensional  potential  vortex,  ccnlrifugal  forces  can  cause  seed 
particles  to  move  radially  away  from  the  votlex  origin.  Con¬ 
sequently,  LV  velocity  measurements  within  vortical  flow- 
nclds^'’  may  be  biased.  For  an  incompressible,  potential  vor- 
lex  flow,  velocity  bias  can  be  computationally  investigated  usina 
a  simplified  equation  of  motion  derived  by  Dring  and  Suo,* 
who  prescribe  particle  position  in  terms  of  three  dimensionless 
parameters  that  are  constant  throughout  the  flow.  For  this  case, 

-  (1) 
ro 

where 


(2) 

(3) 

and 


-  P»Vo</’ 

St  =  — — 

(4) 

IStoM 

For  compressible  flow,  all  three  independent  parameters  vaty 
radially  through  the  vortex.  If  these  varies  are  held  constant 
by  defining  dim  in  terms  of  the  total  density  and  temperature 
of  the  flot^M,  and  trajectories  are  comput^  by  an  equation 
of  motion  which  includes  compressibility  effects,  it  can  be  seen 
in  Fig.  1  that  as  the  initial  Mach  number  of  the  particle  in¬ 
creases  into  the  compressible  flow  regime,  the  particle  trajec¬ 
tory  shows  a  substantial  increase  in  bias.  Consequently,  in  a 
compressible,  vortical  flow,  an  incompressible  analysis  can 
underestimate  the  particle  velocity  bias. 


Approach 

PotcoUal  Vortex 

Although  classical  text  book  derivations  assume  the  poten¬ 
tial  vortex  to  be  an  incompressible  flow  when  deriving  the  ve¬ 
locity  field,’''  the  approach  is  extendible  to  compressible  flow- 
fields.  In  a  compressible  or  incompressible  flow,  the  potential 
vortex  field  is  inotational.  Therefore,  it  is  required  that 


du,  1  /  3ii,\  „ 

U  =  -  -!--(«• - 1  *  0 

>r  r\  »9/ 


<51 


Since  a  potential  vortex  has  no  radial  velocity  component,  and 
since  the  angular  velocity  component  is  only  a  function  of  ra¬ 


preiilble  flow. 


dial  distance,  Eq.  (5)  reduces  to  a  first-order  ordinary  differ¬ 
ential  equation  with  an  integration  constant  relaicd  to  the  vor¬ 
tex  strength  by  the  application  of  Stokes'  theorem: 

r  =  -fu  ■  ds  (6) 


The  resulting  vortex  velocities  are 


u,  = 


2irr 


(7) 


and 


«,  =  0 


(8) 


Since  this  velocity  field  satisfies  the  compreuible  form  of 
the  continuity  equation,  the  flowfleld  is  physically  possible.  If 
die  radial  direction  of  the  vortex  lies  on  the  y-r  plane  of  a 
Cartesian  coordinate  system,  the  two-dimensional  vortical  field 
can  be  superimposed  with  a  constant  normal  velocity  in  the  x 
direction.  The  resulting  three-dimensional  flowfleld  will  then 
have  helical  streamlines  that  rotate  about  the  x  axis.  This  flow- 
field  can  be  used  as  a  simple  model  of  vortices  that  are  shed 
from  aerodynamic  surfaces,  recognizing  that  it  does  not  ac- 
cotmt  for  axial  variations  of  vortex  strength  or  the  development 
of  a  viscous  core. 

Anodier  consideration  for  a  potential  voitex  in  a  compress¬ 
ible  flow  is  the  limit  of  the  mathematical  model.  For  the  in¬ 
compressible  case,  the  temperanire  it  considered  conttaat 
throughout  the  field,  leaving  only  the  infinite  velocity  at  the 
center  of  the  vortex  as  a  limit.  In  a  compiettible  flow,  how¬ 
ever,  the  temperature  decreases  toward  zero  at  the  velocity 
becomes  infinite.  To  examine  the  inner  limit  with  respect  to 
commonly  used  constitutive  equations  for  thermodynamic  clo¬ 
ture.  it  it  best  to  temporarily  decouple  the  constam  axial  ve¬ 
locity  from  the  vortical  flow  plane.  Since  the  flow  can  be  con¬ 
sidered  isentropic  and  adiabatic,  the  axial  velocity  can  be  used 
to  relate  the  total  temperature,  pressure,  and  density  of  the 
flowfleld  to  the  static  values  in  the  outer  vortex  region  which 
approaches  zero  velocity  in  the  vortical  flow  plane.  These 
'outer*  values  can  then  be  considered  as  stagnation  values  for 
the  decoupled  two-dimensional  vortex  flowfleld.  Within  the 
region  of  the  flowfleld  where  the  perfect  gas  law  and  calori- 
cally  perfect  gas  assumption  hold. 


and 


where 


E.  = 

P. 


1 


r 

8»V’c,r, 


(9) 

(10) 


(II) 


(12) 


The  thermodynamic  profiles  in  terms  of  the  transformed  ra¬ 
dial  coordinate  are  shown  in  Fig.  2.  For  air,  the  stagnarkm 
temperature  within  the  vortical  plane  can  be  as  high  as  470  K, 
bas^  on  an  arbitrary  maximum  of  2%  error  in  the  calorically 
perfect  gas  assumption.  The  inner  limit  of  the  nwdel,  based 
on  a  maximum  3%  error  in  Sutherland’s  viscosity  law,  it  at 
radial  distance  from  the  vortex  origin  where  the  ttagiwtion 
temperature  has  decreased  to  167  K. 

The  corresponding  range  of  P,  extends  fiom  near  vacuum 
to  an  upper  limit  defined  by  the  validity  of  the  perfect  gat  law. 
For  a  maximum  2%  error,  any  F,  up  to  11.2  atm  is  valid  over 
the  entire  range  of  T,. 
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Fit-  1  Tbcrmodynamk  proHIci  through  a  poUnllai  rortcx  In  a 
conprcniblc  flow. 


Partkk  £<|«atioa  at  Modon 

The  equation  set  used  to  describe  the  particle  motion  is  de¬ 
veloped  by  Nicholt,*  and  includes  the  influences  of  Stokes* 
drag,  resultant  pressure,  apparent  mass,  gravitational  forces, 
and  the  Basset  history  integral.  For  application  to  the  vortical 
flow,  the  gravitational  term  and  the  Basset  history  integral  have 
been  excluded.  The  Basset  term  can  be  significant  when  the 
particle  density  is  of  the  same  order  of  magnitude  as  the  fluid 
density,  or  when  the  particle  relative  acceleration  is  large  com¬ 
pared  with  the  particle  relative  velocity.  In  this  case,  as  in 
most  numerical  solutions,  the  Basset  term  is  neglected  because 
the  numerically  evaluated  integral  is  found  to  be  either  neg¬ 
ligible,  highly  unstable,  or  potentially  biased  by  the  accumu¬ 
lation  of  roundoff  and  truncation  errors.  The  resulting  differ¬ 
ential  equation  for  a  general  steady-state  flow  can  be  wrinen 
as 


dv,  -3C^ 
dr 


4rf  (m,  +  nv/2) 
m, 


■  W,  H>  + 


m, 


(m,  +  mf/2) 


(m,  +  mf/2)  '  dx, 

-{u.  +  ^vv} 

dr  I  40  J 


(13) 


where 


tv,  =  V,  -  u,  - 


(14) 


and  each  of  the  fluid  properties  and  derivatives  are  evaluated 
with  respect  to  the  particle  location. 

The  drag  coefTicienl  is  evaluated  by  the  drag  law  of  Crowe”®. 

Cd  =  (Cd,  -  2)  exp|-3.07r''’(Af,/lfe,)C(Ke,)) 


Crowe's  drag  law  is  valid  over  a  wide  range  of  relative 
Reynolds  and  Mach  numbers,  and  includes  the  effects  of  in¬ 
ertia,  compressibility,  rarefaction,  and  heat  transfer. 

Nichols*  shows  Ihitt  the  particle  equation  of  motion  becomes 
stiff  for  small  particle  diameters.  Consequently,  his  soluiioru 
for  particle  behavior  through  shock  waves  ate  piecewise,  using 
Adams'  technique  where  possible,  and  restarting  the  solution 
with  slower  Runge-Kutta  integration  when  local  stiffnew  causes 
Adams’  method  to  diverge.  To  avoid  the  necessity  of  piece- 
wise  solutions  while  maintaining  computational  efficiency,  so¬ 
lutions  presented  here  are  obtained  by  4th  order  Runge-Kutta- 
Fehlberg  (RKF)  numerical  integration.  The  RKF  method  is 
ideally  suited  for  this  situation  since  it  estimates  the  local  stiff¬ 
ness  of  the  equation  at  each  time  step,  aiul  adjusts  the  step  size 
accordingly.  Results  include  a  gen^  consideration  of  trajec¬ 
tories  and  velocity  bias  within  the  potential  vortex  naodel,  spe¬ 
cific  ai^lication  of  die  model  to  an  experimental  LV  survey, 
and  the  use  of  a  compuutional  fluid  dynamics  (CFD)-gener- 
ated  vortical  flowfield  in  place  of  the  model. 

Results 

Dimeiuloiial  Analysis 

The  particle  equation  of  motion  described  by  Eq.  (13)  re¬ 
quires  the  coupled  solution  of  three  first-order  differential 
equations  initialized  by  particle  position  and  velocity  widt>  a 
specified  vortical  field.  To  reduce  the  equation  set  for  analysis, 
consider  the  slightly  more  specific  case  of  an  LV  seed  particle 
traveling  downstream  in  ai^w  (in  the  x  direction),  and  at 
time  zero,  it  becomes  surrounded  by  a  potential  vortex.  If  the 
particle  it  traveling  with  the  tame  initial  constant  downstream 
velocity  as  die  flow,  the  equation  of  motion  for  the  x  direction 
decouples  from  the  vortical  flow  in  the  y-z  plane.  Since  there 
are  no  initial  forces  on  the  particle  in  the  x  directioo,  the  equa¬ 
tions  show  that  the  downstream  particle  velocity  will  remain 
constant.  Additiorully,  since  the  vortex  is  axisymmetric,  one 
initial  condition  parameter  can  be  eliminated  hy  choosing  yo 
as  the  initial  particle  position,  and  setting  zo  to  mid.  Also, 
since  the  particle  is  assumed  to  have  no  initial  velocity  other 
than  in  the  X  direction,  both  of  the  initial  velocity  parameters 
in  the  vortical  plane  can  be  set  to  zero. 

Analysis  of  the  reduced  equation  set  in  the  vortical  plane 
coupled  with  the  various  equations  of  closuie  require  the  fol¬ 
lowing  dimensional  inputs: 


“0- P- r„  5,  </.  p,,  yo  (19) 

where  5  is  Sutherland's  constant  for  air,  and  the  outer  viscosity 
can  be  evaluated  by  Sutherland’s  law  for  the  input  value  of 
T,.  These  dimensional  inputs  describe  the  dependent  variables; 

y.  r,  V,.  V,.  r  (20) 

Application  of  the  Buckingham  11  Theorem'  *  reduces  the 
dimensional  equation  set  to  a  function  in  terms  of  five  dimen¬ 
sionless.  independent  parameters: 

T,  ^  p.ihd  ^  Ya 
S  HT;  p.  ’p/d 


+{«Af,)/(y’'*W»)}  exp{  -fle,/(2A#,)}  +  2  (15) 

where 

Cd,  =-^(1  +  O.lSWfri'’)  (16) 

Ke, 

log,oC{/fe,)  =  l.25{l  +  ianh(0.77  log.o^e,  '  1.92)}  (17) 

and 

/!(«,)  =  {2.3  i.7(r,/r/)''’) 

-  2.3  tanh(  1,17  log, oM»)  (1*) 


Additional  trial-and-error  analysis  shows  that  the  effect  of 
the  parameter  involving  Sutherland's  constant  is  negligible,  and 
that  the  density  ratio  and  initial  distance  parametets  can  be 
combined.  The  resulting  dimensionless  equation  set  is 


L  L  ^  ^  T=/(— 

yo  yo  “o’  “o'  \kT.'  p,  ’dp,] 


(22) 


where 


T  =  ruo/yo  '23) 

Two  of  the  dependent  dimensionless  parameters  in  Eq.  (22) 
describe  the  particle  velocities  in  the  vortical  plane.  However, 
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describe  the  particle  velocities  in  the  vortical  plane.  However, 
for  application  to  LV,  it  is  the  bias  in  these  quantities  that  is 
important.  Therefore,  define 


V,  -  He 

%£RR,  =  - - 7-775  X 

100 

(24) 

V.  -  u. 

%EM,  =  X 

(u’  +  nf)'-’ 

100 

(25) 

Since  these  two  error  parameters  can  be  written  as  functions 
of  the  dimensionless  groups,  they  can  be  used  as  dependent 
dimensionless  parameters  themselves. 

The  three  independent  dimensionless  groups  can  be  fully  de¬ 
fined  by  inputs  of  outer  temperature  and  density,  initial  flow- 
field  velocity  and  particle  position,  and  the  particle  density  and 
diameter.  In  terms  of  dimensionless  variables,  these  groups 
can  be  defined  as 


tr,  = 

(26) 

iTi  =  Potkiditt, 

(27) 

wj  =  (yoP.)/{dp,) 

(28) 

The  six  independent  dimensional  parameters  and  the  vis¬ 
cosity  from  Sutherland's  law  are  bounded  for  LV  purposes  in 
Table  I.  The  corresponding  ranges  of  the  three  »  parameters 
are  also  shown. 

Particle  trajectories  in  the  y-z  plane  for  a  single  combination 
of  tr  groups  are  shown  in  Fig.  3.  For  this  case,  the  particle 
has  moved  to  more  than  twice  its  initial  radial  position  after 
three  helical  revolutions.  Although  the  trajectory  alone  dbes 
not  quantify  the  velocity  bias,  it  does  show  the  difficulty  in 
maintaining  seeded  flow  near  the  vortex  origin. 

For  the  same  combination  of  ir  groups,  the  measurement 
error  that  will  occur  in  the  y  and  z  velocity  components  as  a 
function  of  dimensionless  time  is  shown  in  Fig.  4.  In  general, 
the  magnitude  of  velocity  bias  decreases  from  100%  at  the 
initial  condition  toward  0%  as  the  particle  moves  outward  from 
the  vortex  origin  to  a  radius  where  it  accurately  follows  the 
fluid  streamlines.  Since  the  velocity  normal  to  the  vortical  plane 
is  constant,  values  of  r  correspond  linearly  to  particle  trans¬ 
lation  along  the  x  axis. 

For  the  ranges  of  parameters  selected  in  Table  I ,  the  radial 
particle  position  and  corresponding  time  at  which  the  total 
magnitude  of  velocity  bias  decrease:^  to  3%  is  shown  in  Figs. 
3a- 5c.  On  the  log-log  scale,  the  curves  of  radial  position  are 
shown  to  decrease  nearly  linearly  towards  conditions  where 
any  radial  displacement  of  the  particle  is  negligible.  In  terms 
of  the  dimensional  parameters,  this  decrease  corresponds  to 
smaller  particle  diameters,  an  initial  position  further  from  the 
vortex  origin,  or  lighter  particles  relative  to  the  fluid  density. 

The  corresponding  curves  for  t  also  decrease  as  ir,  in¬ 
creases,  and  include  a  sharp  drop  at  the  point  where  radial 


translation  becomes  negligible  In  general,  damping  the  ve¬ 
locity  bias  from  the  initial  conditions  is  much  faster  than  any 
required  radial  translation  of  the  particle.  Consequently,  the 
portion  of  the  curve  above  the  dip  is  fur  particles  that  must  be 
translated  further  away  from  the  vortex  origin,  while  the  por¬ 
tion  of  the  curve  below  the  dip  conesponds  to  particles  that 
are  brought  from  rest  to  the  fluid  velocity  at  the  initial  radial 
particle  position. 

Application  of  the  Potential  Vortex 

As  a  specific  application,  the  potential  vortex  is  used  to  model 
an  LV  measured  supersonic  vortex  shed  from  a  7S-deg  delta 
wing  at  20-deg  angle  of  attack.  The  two-dimensional  velocity 
vector  field  at  80%  chord  is  shown  in  Fig.  6,  and  a  compre¬ 
hensive  description  of  the  experimental  test  and  results  can  be 
found  in  Ref.  3.  From  the  known  stagnation  conditions  and 
fteestream  Mach  number,  the  outer  vortex  temperature,  den¬ 
sity,  and  axial  velocity  are  determined,  and  the  density  of  the 
se^  corresponds  to  the  10  cS  silicon  oil  that  was  used  for  the 
test. 

Three-dimensional  LV  measurements  were  made  in  dte  plane 
perpendicular  to  the  mrxlel  surface.  However,  it  was  found 
that  by  applying  a  rotational  transformation  to  the  data,  lines 
emanating  normal  trom  the  individual  two-dimensional  veloc¬ 
ity  vectors  tend  to  converge  towards  a  point.  This  point  is  es¬ 
timated  to  be  the  vortex  origin.  The  transformed  velocity  field 
shown  in  Fig.  6  conesponrls  to  a  ‘vortex  angle  of  attack*  of 
8  deg  from  the  fteestream  flow. 

The  erxtesponding  values  of  vortex  strength  decrease  from 
the  outer  edge  of  the  wing  towards  the  delta  wing  centerline. 
For  this  example,  a  nominal  weighted  average  value  of  -30 
m'/s  is  used. 

Using  the  inrrer  limit  of  the  potential  vortex  model  as  an 
initial  position,  the  corresponding  particle  trajectories  for  var- 
irxis  sizes  of  seed  ate  shown  in  Fig.  7.  Because  of  the  high 
speed  of  the  axial  flow,  the  particles  only  travel  abrxit  rxie- 
third  of  a  helical  revolution  from  the  nose  of  the  wing  to  the 
80%  chord  location.  However,  the  difference  in  individual  tra¬ 
jectories  it  significant,  showing  a  nearly  circular  path  for  the 
O.l-srm  particle  while  the  5.0-fzm  particle  increases  its  radial 
disUnce  from  the  vortex  origin  by  mote  than  40%. 

The  radial  position  of  particles  at  the  80%  chord  location  as 
a  continuous  fiinction  of  seed  diameter  is  given  in  Fig.  8.  Over 
the  included  range  of  0. 1-3.0  ftm  particles,  the  curve  shows 
that  the  increase  in  radial  translation  is  fairly  linear. 

The  corresponding  curve  of  the  velocity  bias  at  this  chord 
location  is  given  in  Fig.  9.  This  curve  is  also  fairly  linear, 
showing  a  steady  increase  in  bias  from  less  than  0.3%  for  0. 1- 
run  particles  to  nearly  28%  for  S.O-rrm  seed.  This  curve  shows 
that  seed  diameters  of  the  polydispersed  silicon  oil  must  be 
less  than  0.833  fon  in  order  to  maintain  a  velocity  bias  of  less 
than  3.0%. 

Finally,  the  measurement  error  for  three  candidate  particle 
sizes  as  a  function  of  chord  position  is  shown  in  Fig.  10.  In 
each  case,  the  particle  bias  at  zero  percent  chord  is  100%, 


Table  I  Chosen  ranges  of  dimensional  and  corregrondlng  dlracnsloolcis  parameters 


Parameter 

Range 

Reasoning 

r„  K 

176-470 

The  lower  value  it  aibitrahly  chosen  so  ihai  M  2  0.5  at  the  inner  limit  of  the  vortex  model.  The 
upper  value  is  based  on  a  maximum  2%  error  for  the  caiorKally  perfect  gas  assumption. 

P,.  kg/m’ 

0.0375-4  011 

Dependent  on  T..  this  range  corresponds  to  the  minimal  range  of  0.05  atm  <  Po  <  2  atm.  which 
is  an  arbitrary  range  of  interest. 

kg/ms 

1  I9E-5-2  56E-5 

From  Sutherland's  viscosity  law  for  the  range  of  To. 

Ud.  m/$ 

1(7-780 

The  lower  value  is  arbitrary,  the  hi^  value,  dependent  on  To,  is  the  velocity  at  the  inner  limit  of 
the  vonex  model. 

P,.  kg/m’ 

260-3960 

This  range  covers  a  wide  variety  of  LV  seed,  from  hollow  glass  spheres  to  AljO). 

d.  ftm 

0  2-100 

This  is  a  typical  range  for  LV  seed  diameters. 

>0.  m 

0.001-2  0 

This  is  a  probable  range  of  interest  for  wind-tunnel  applications 

^1 

0.0273-2  144 

rr,.  try.  uid  n,  are  based  on  appropnaie  combinations  of  dimensional  parameters. 

0.00293-12221 

»i 

9  47  X  10' ’-1 54269 
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Fig.  4  Rcprticntadvc  (line  hWorr  of  rdoclly  Mai  for  one  pm- 
liHt  n  gnwp. 


since  (he  particles  have  no  initial  velocity  in  the  vortical  plane. 
By  the  time  the  particles  move  downstream  only  2-5%  chord, 
however,  the  initial  bias  is  completely  damped.  The  remaining 
bias  is  due  to  the  tangential  momentum  of  the  particle  in  the 
vortical  Tield,  and  slowly  decays  as  (be  radial  distance  of  the 
seed  from  the  vortex  origin  increases. 

Consequently,  this  example  predicts  that  for  seed  particles 
of  leu  than  0.833  pm  in  diameter,  reasonably  unbiased  LV 
measurements  can  be  made  within  the  radial  limits  of  the  po- 
temial  vortex  model  at  any  axial  station  downstream  of  about 
10%  chord.  Since  the  inner  bound  of  the  model  happens  to  be 
the  approximate  bound  at  which  LV  measuremenu  were  ac¬ 
tually  taken,  the  unmodeled  viscous  inner  core  of  the  vortex 
is  inconsequential  to  this  analysis.  Furthermore,  it  wu  found 
that  any  possible  initial  velocities  applied  to  (he  seed  particles 
were  quickly  damped,  making  (he  results  nearly  independent 
of  the  simplified  set  of  initial  conditions  used. 

AppOcallan  of  a  Compalallonai  FhiM  Dynamlcs-Ccncratcd 
Vortical  FMd 

As  a  second  qiplkation  of  the  particle  equations  of  motion, 
the  potential  vortex  model  is  replaced  with  a  three-dimensional 
Navier-Stokes  solution  of  a  delta  wing  flowfield.  For  (he  ume 
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model  geometry  at  slightly  different  flight  conditions  dian  the 
experimental  case,  two-dimensional  crou  sections  of  a  portion 
of  (he  delta  wing  flowfield  calculated  by  Webster  and  Shang" 
are  shown  in  Fig.  1 1 .  As  represented  by  the  IS  and  40%  chord 
locations,  the  general  form  of  the  primary  and  secondary  vor¬ 
tices  appear  to  vary  only  slightly  with  chord,  although  (he  size 
of  the  structures  varies  proportionally  with  the  wing  span. 

As  an  atten^  to  seed  the  inner  core  of  the  vortex,  seed 
pmticles  were  initially  positioned  in  the  vicinity  of  the  leading 
edge  of  the  model.  It  was  found,  however,  that  the  vertical 
momemum  of  the  particles  from  the  freestream  angle  of  attack 
is  too  great.  The  particles  enter  the  lower  poition  of  the  form¬ 
ing  vortex,  and  continue  to  move  vertic^ly  upward  through 
the  cote  and  into  the  outer  portion  of  the  vortex  before  the 
initial  vertical  momentum  is  dissipated.  By  trial-and-enor,  it 
was  found  that  the  seed  which  peneuates  deepest  into  the  vor¬ 
tex  cote  are  particles  that  are  swept  around  the  edge  of  the 
wing  from  (he  underside. 

Trajectories  of  seed  particles  that  are  initially  positioned  along 
the  edge  of  the  wing  at  15%  chord  are  shown  in  Fig.  12.  As 
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Fig.  12  Two-dliatMiaaal  pnjcctiaa  «f  rilkaa  oH  Ktd  partkk 
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the  particles  travel  downstream,  they  are  swept  around  the  pri¬ 
mary  vortex.  The  smallest  particle  shown.  0. 1  iim.  completes 
an  entire  helical  revolution  prior  to  reaching  the  40%  chord 
station,  and  maintains  a  velocity  bias  of  less  than  3%.  The 
2.0-(un  particle,  on  the  other  hand,  is  centrifuged  out  of  (be 
vortex,  and  has  a  final  velocity  bias  of  over  30%. 

Trajectories  for  0.2-srm  particles  at  different  initial  chord 
locations  along  the  edge  of  the  wing  are  shown  in  Fig.  13.  In 
each  case,  the  partKies  are  (raveling  helically  around  the  i/,-; 
mary  vortex,  which  is  growing  chordwise.  The  particle  (hat 
starts  at  the  S%  chord  location  is  shown  to  travel  otkc  around 
the  primary  vortex,  and  then  it  breaks  off  to  circle  the  sec¬ 
ondary  vortex.  As  that  particle  moves  beyond  the  40%  chord 
location  shown,  it  passes  back  from  the  secondary  to  the  pri¬ 
mary  vortex. 

If  the  40%  chord  positions  of  each  of  the  particles  shown 
in  Fig.  13  are  connected,  the  resulting  curve  will  be  the  inner 
core  bound  for  the  seed  found  at  (hat  profile  location.  This 
curve  is  shown  in  Fig.  14  for  three  particle  diameters.  As  would 
be  expected,  the  size  of  the  unseeded  inner  core  grows  with 
the  size  of  the  seed.  The  0.5-/tm  particles  travel  slightly  less 
than  one  helical  revolution,  and  have  an  average  velocity  bias 
at  40%  chord  of  about  10%.  The  0.2-/tm  particles  have  an 
average  velocity  bias  of  around  S%,  and  the  seed  (hat  starts 
near  the  leading  edge  is  shown  to  be  circling  the  secondary 
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vortex  at  the  40%  chord  location.  The  0  l-ftm  seed,  with  an 
average  velocity  bias  of  about  3%,  is  entirely  within  the  pri¬ 
mary  vortex  at  the  40%  chord  measurement  station.  Although 
the  specific  submicron  particle  size  distribution  of  the  experi¬ 
mental  study  is  nut  known,  the  size  and  shape  of  the  unseeded 
core,  as  determined  by  the  LV  surveys  and  laser  light  sheet 
visualization,  correspond  qualitatively  with  the  analytical 
ptediction. 


ConduskNU 

Three  contributions  were  made  by  this  study.  First,  a  three- 
dimensional  potential  vortex  flowTield  was  developed  and 
bounded  for  use  as  a  mathematical  model.  Second,  curves  were 
presented  that  can  be  used  to  aid  in  the  design  and  analysis  of 
LV  measurements  in  vortical  flows  over  a  large  practical  range 
of  interest.  Third,  (he  particle  equations  of  motion  were  cou¬ 
pled  with  a  CFD-dcriv^  fiowfieM,  showing  (he  capability  to 
numerically  track  the  position  and  velocity  bias  of  LV  seed 
through  highly  complex  flows. 
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Allhough  the  polenlial  i>  a  very  simple  model  of  Ihe 

della  wing  vociex  struei'  .e  it  gave  valuable  insight  into  Ihe 
behavior  of  the  silic.  >il  seed  in  the  wind-tunnel  environ¬ 
ment  It  eslimaieJ  'he  required  panicle  diameter  to  within  Ihe 
same  order  of  rr..  ^.ulude  as  using  the  CFD  solution,  and  offers 
results  for  the  entire  class  of  solutions  in  terms  of  three  di- 
inensionless  variables. 

Tracking  Ihe  panicles  through  the  CFD-generated  flowTield 
offers  the  greatest  potential  to  investigate  seed  density  distri¬ 
butions  and  panicle  bias  in  any  LV  study  which  is  coupled 
with  a  CFD  effon.  Using  the  RKF  melh^  of  numerical  in¬ 
tegration  to  automatically  adjust  the  lime  step  based  on  eq'sa- 
lion  stiffness  makes  Ihe  method  numerically  stable  and  efliciesl. 

Both  methods  predict  that  the  velocity  bias  within  the  delta 
wing  vones  increases  rapidly  with  panicle  diameter.  Conse¬ 
quently.  panicle  size  must  be  carefully  controlled  in  order  to 
minimize  panicle  dynamic  bias.  Future  effons  will  use  Ihe  po¬ 
tential  vonex  model  to  design  expenmenis  and  analyze  data 
for  vonical  flows  that  do  not  have  a  companion  CFD  solution. 
When  a  CFD  solution  exists,  it  will  be  used  with  the  panicle 
equation  of  motion  as  a  more  precise  method  of  quantifying 
the  dynamic  measurement  bias. 
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APPENDIX  C:  COMPUTATIONAL  FLUID  DYNAMICS  RESULTS 


This  appendix  contains  details  of  the  computational  fluid  dynamics  study  of  the 
supersonic  vortex  flowfield.  It  is  a  reproduction  of  an  AIAA  Journal  article  (Ref  13)  written  as 
a  work  of  the  Aeromechanics  Division,  Flight  Dynamics  Directorate,  Wright  Laboratory, 
Wright-Patterson  Air  Force  Base. 
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Thin-Layer  Full  Navler-Stokes  Simulations 
over  a  Supersonic  Delta  Wing 
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Steady  flowflekb  dcaeriblaf  the  flow  over  a  7S-drt  tarcpi  delta  wiat  at  -  \.9S  aad  4.41  x  II*  were 
•imulatcd  at  a  -  M  aad  M  dej  aiing  two  eoaipaler  codea.  CoBipailsoat  were  aude  betweea  cakulatloat  aaiag  1) 
the  lamlaar  aad  tarbitleat  Ibla-layer  Navler-Stokee  cqaadoat,  aad  2)  teailnar  tbla-layer  aad  Navler-Stokea  cqaa- 
tioas.  At  a  •  26  def,  eacb  eqaatkia  let  captarei  tbe  eiaeallal  Mnietare  of  the  flow  aad  the  differeacei  betweea 
their  resalu  are  nlaor.  Naaierical  ramiti  geaerated  by  a  grid  reflaeawat  itudy  eiblbltcd  oaly  a  ailaor  iBprove- 
aieat.  The  Navler-Stokei  egaatloai  were  aied  to  calealate  tbe  flow  at  a  >  30  deg.  A  reglea  of  renned  flow  aloag 
tbe  lurface  aear  tbe  traHlag  edge  wet  obtcrred.  Tbe  eeetkal  exteat  of  Ibit  regioa  was  aiach  aatflcr,  bat  tbe 
apilnaai  propagatloa  exteaded  aiach  farther  tbaa  the  thladayce  daialatloa. 


Nomenciafure 

Cp  >coerricient  of  pressure 

L  <=  root  chord  length 

M  =  Mach  number 

Re  s  Reynolds  number 

U,  K  W  =  velocity  components  in  coordinate  directions 
jf.y.z  s  coordinate  directions 

Y  =■  maximum  span  at  a  given  x  location 

y'  s  nondimensional  distance  from  surface, 

rv/rp/p/e) 

a  »  angle  of  attack 

d  ••  leading-edge  bevel 

r  =  sweep  angle 

p  « density 

r  <=  surface  shear  stress 

oD  3  freestream  value 

0  surface  value 

Inirotluctitin 

HE  flow  about  modem  high-speed  aircraft  with  a  delta 
planform  at  angle  of  attack  is  characterized  by  the  pret¬ 
ence  of  large  spiraling  vortices  on  the  lee  side  of  the  wing. 
In  many  cases,  these  vortices  are  the  primary  structure  in  the 
flow  affecting  the  performance  of  the  aircraft.  These  vtKticcs 
readily  form  at  low  angles  of  attack.  With  increasing  angle  of 
attack,  secondary  and  tertiary  vortices  form  due  to  boundary- 
layer  separation.  Flows  over  slender  delta  wings  have  been 
studied  experimentally  by  Monnerie  and  Werle,',  Hummel.’ 
Miller  and  Wood,’  Stallings  and  Lamb,*  and  Guyton,’  among 
others;  they  have  been  studied  numerically  by  Rizzetta  and 
Shang,*  Buter  and  Rizzetta,’  Thomas  and  Newsome,'  to 
name  but  a  few.  There  have  been  numerous  investigations  of 
other  delta-based  planforms.  These  include,  but  are  not 
limited  to,  delta  wings  with  canard,*  crank  delta  wings," 
pitching  delta  winp,"  and  rocking  delta  winp.” 

In  a  previous  paper,”  the  authors  presented  numerical  sim¬ 
ulations  of  vortical  flow  over  a  delta  wing  (T  •  7}  deg)  at  four 
angles  of  attack  (10,  20,  30,  and  33  deg)  at  a  freestream  Mach 
number  of  1.93  and  a  Reynolds  number  of  4.4Sx  10*.  The 
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governing  equations  used  were  the  unsteady,  compressible, 
three-dimensional,  thin-layer  Navier-Stokes  equations  written 
in  Reynolds  averaged  form.  The  computer  code  used  to  solve 
this  set  of  equations  was  ARC3D,  widch  is  fully  described  by 
Pulliam.”  This  code  uses  the  a^oximately  factorized  im¬ 
plicit  Beam-Warming  algorithm  and  is  diagonalized  to  acce¬ 
lerate  convergence. 

The  previous  results  captured  the  essential  structure  of  the 
vortical  flow  on  the  Ice  side  of  the  delta  wing.  At  a  -  10  and  20 
deg,  the  results  compared  very  favorably  with  experimental 
results.  At  aw  30  and  33  deg,  a  large  region  of  reversed  flow 
was  observed  near  the  trailing  t^e.  This  reversed  flow  was 
contained  in  the  low-pressure  region  associated  with  the  secon¬ 
dary  vortex.  It  extended  from  the  trailing  edp  forward  to 
679b  of  the  root  chord  for  a  •  30  deg  and  forward  to  36%  of 
the  root  chord  at  a  »  33  deg. 

This  work  left  three  areas  in  which  the  results  were  not 
satisfactory.  The  first  problem  area  was  the  assumption  of 
laminar  flow.  The  high  Reynolds  number  (4.46  x  10*  based  on 
root  chord  length)  would  indicate  the  flow  to  be  turbulent. 
However,  to  focus  the  range  of  issues  in  the  previous  work, 
turbulence  was  not  modeled.  The  second  issue  was  the  use  of 
the  thin-layer  assumption.  The  thin-layer  assumption  is  con¬ 
sidered  appropriate  for  vortical  flows;  however,  its  utility  is 
unknown  in  the  reversed  flow  regioa  where  stress  components 
in  tbe  spanwise  direction  may  not  be  negligible.  The  final 
problem  area  was  tbe  need  for  a  grid  resolution  study  that  was 
not  performed  in  the  previous  work. 

Addressing  eacb  of  these  issues  is  the  objective  of  this 
paper.  In  order  to  do  this,  a  second  computer  code  was  used  to 
solve  the  Navler-Stokes  equations.  This  code  (FDL3D-I),  writ¬ 
ten  by  Visbal,”  it  based,  like  ARC3D,  on  the  approximately 
factorized  implicit  Beam-Warming  algorithm;  however,  it  has 
several  differences  from  ARC3D.  First,  the  code  can  solve 
either  the  full  or  thin-layer  Navier-Stokes  equations.  Second, 
the  code  it  not  diagonalized  but  retains  the  block  tridiagonal 
form;  thus,  it  it  first-order  accurate  in  time.  In  addition  to 
these  major  differences,  there  are  several  smaller  differences 
between  the  two  codes.  Since  FDL3D-I  retains  the  block  tri¬ 
diagonal  structure,  the  implicit  damping  terms  are  second 
order.  BecauK  of  its  scalar  pentadiagonal  form,  ARC3D  has 
both  second-  and  fourth-order  terms.  In  addition,  the  metric 
terms  at  boundaries  in  ARC3D  are  first-order  accurate, 
whereu  in  FDL3D-I,  they  are  second  order. 

In  order  to  compare  the  codes  as  carefully  as  pouible,  the 
damping  terms  were  maintained  as  nearly  identical  as  possible. 
The  coefficients  for  the  second-  and  fourth-order  explicit  and 
second-order  implicit  were  the  same.  To  match  the  fourth- 
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order  implicit  term  from  AR03D,  additional  second-order  im¬ 
plicit  damping  was  added  to  FDL3D.  This  results  in  the  codes 
solving  slightly  different  sets  of  equations,  and  small  differ¬ 
ences  in  the  results  were  to  be  expected.  One  method  to  avoid 
this  difference  would  have  been  not  to  use  fourth-order  damp¬ 
ing;  however,  the  solution  became  unstable. 

The  resolution  of  the  issues  just  outlined  was  accomplished 
in  the  following  manner.  First,  the  effects  of  turbulence  were 
evaluated  using  ARC3D.  This  was  done  since  ARC3D  has  an 
implementation  of  the  Baldwin-Lomax  turbulence  model, 
which  is  currently  unavailable  in  FDL3D-I.  Second,  FDL30-I 
was  used  to  calculate  the  flowfield  about  a  delta  wing  using  the 
thin-layer  approximation  to  validate  FDL3D-I.  This  calcula¬ 
tion  us^  the  same  conditions  and  grid  as  the  previous  work  at 
a  =  20  deg.  These  results  were  compared  to  the  previous  results 
as  well  as  to  experimental  data.  Third,  the  full  Navier-Stokes 
equations  were  solved  using  FDL3D-I  and  compared  to  the 
thin-layer  results,  such  that  observed  variations  in  the  results 
could  be  attributed  to  the  use  of  thin-layer  or  full  Navier- 
Stokes  equations  and  not  to  differences  between  the  two 
codes.  Fourth,  the  effects  of  grid  resolution  were  investigated 
using  the  Navier-Stokes  equations.  The  spanwise  resolution 
over  the  delta  wing  and  the  normal  resolution  near  the  upper 
surface  were  increased  by  a  factor  of  2. 

All  of  this  work  was  done  at  a  =  20  deg,  a  condition  that  did 
not  display  significant  reversed  flow  in  the  previous  work.  In 
order  to  determine  the  effects  of  using  the  full  Navier-Stokes 
equations  on  the  reversed  flow  region,  the  a  -  30-deg  case  was 
solved  using  FDL3D-I.  This  angle  of  attack  produced  a  signi¬ 
ficant  region  of  reversed  flow  in  the  previous  results. 

Analysif 

The  vortical  flow  about  a  delta  wing  (r-73  deg)  was 
simulated  numerically  at  a  =  20  deg  and  30  deg,  Af.  =  1 .95 
and  4.48  X  10*.  The  delta  wing  configuration  had  a  length  of 
33.9  cm,  a  thickness  of  1 .9  cm,  and  a  leading-edge  bevel  0  =  }S 
deg.  Figure  I  is  a  sketch  of  the  delta  wing  with  relevant  param¬ 
eters  labeled.  This  configuration  was  chosen  to  duplicate  a 
parallel  experimental  study  that  is  underway  at  the  present 
time. 


Fig.  I  Sketrh  of  Mta  wtag. 


The  governing  equations  were  taken  to  be  the  unsteady, 
compressible,  three-dimensional  Navier-Stokes  equations. 
The  perfect  gas  law  was  used  to  relate  the  pressure,  tempera¬ 
ture  and  density,  whereas  Sutherland’s  equation  was  us^  to 
calculate  the  molecular  viscosity.  Two  computer  codes  were 
used  to  solve  this  set  of  equations.  The  first,  ARC3D,  uses  the 
thin-layer  approximation.  The  second,  FD1.3D-I,  has  the  op¬ 
tion  of  solving  either  the  full  or  thin-layer  Navier-Stokes  equa- 


Flg.  1  OalHat  ol  Ike  grM  ttraetart:  i>  global  •Irarlort;  b)  oolboard 
tig;  c)  iraRlag  edge. 
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tions.  The  local  time  step  for  each  code  was  scaled  by 
(I.O-t-^)'',  where  J  is  the  Jacobian  of  the  coordinate 
transformation. 

The  grids  for  these  calculations  were  generated  alge¬ 
braically.  The  grid  may  be  described  as  having  an  H-H  struc¬ 
ture.  This  type  of  grid  was  chosen  to  eliminate  singular  metrics 
at  the  leading  edge  and  to  facilitate  modeling  of  the  blunt  trail¬ 
ing  edge.  Figures  2  show  several  planes  that  outline  the  global 
grid  structure  (Fig.  2a).  The  planes  are  the  following:  the  sym¬ 
metry  plane  (labeled  A  in  the  figure),  the  plane  that  contains 
the  upper  surface  of  the  delta  wing  ^beled  B),  and  the  plane 
that  is  located  at  the  downstream  edge  of  the  domain  (labeled 
C).  In  addition,  Figs.  2  show  close-ups  of  the  outboard  tip 
(Fig.  2b)  and  the  trailing  edge  (Fig.  2c),  detailing  the  H-H 
structure.  The  trailing  edge  was  rounded  in  order  to  smooth 
the  grid.  The  normal  grid  spacing  at  the  surface  increased  pro¬ 
portionally  with  the  distance  from  the  apex  to  a  maximum  of 
10  ’  X  Z.  at  the  trailing  edge.  The  grid  was  clustered  in  the 
streamwise  direction  at  the  leading  and  trailing  edge  with  the 
spacing  on  the  surface  varying  by  almost  three  orders  of  mag¬ 
nitude.  A  more  detailed  discussion  of  the  grid  structure  was  re¬ 
ported  previously.’’ 

The  coordinate  system  had  its  origin  at  the  apex  with  the  x 
axis  chordwise  along  the  upper  surface  in  the  plue  of  symme¬ 
try.  They  axis  was  normal  to  the  plane  of  symmetry  and  in  the 
spanwise  direction.  The  z  axis  was  normal  to  the  upper  surface 
and  in  the  plane  of  symmetry. 

The  primary  grid  had  80  grid  points  in  the  x  direction  with  S 
ahead  of  the  apex  and  IS  in  the  wake.  In  the y  direction,  there 
were  49  grid  points  inboard  of  the  leading  edge  and  31  out¬ 
board  for  a  total  of  80  grid  points.  In  the  z  direction,  there 
were  80  grid  points  with  32  below  and  48  above  the  delta  wing. 
This  resulted  in  a  grid  with  a  total  of  312,000  grid  points.  This 
grid  was  used  for  the  comparisons  between  ARC3D  and 
FDL3D-I,  the  comparisons  between  the  full  and  thin-layer 
Navier-Stokes  equations  using  FDL3D-I,  and  for  the  turbu¬ 
lent  calculation  with  ARC30. 

For  the  grid  refinement  study,  a  second  system  had  an  addi¬ 
tional  49  grid  points  in  the  spanwiK  direction.  A  new  grid 
point  was  inserted  between  each  grid  point  on  the  surface.  In 
the  z  direction,  a  new  point  was  inserted  between  each  of  the 
first  21  points  above  the  upper  surface.  This  resulted  in  a  grid 
with  1 ,032,(XI0  grid  points. 

The  final  grid  was  based  on  the  primary  mesh  system  with 
an  additional  30  points  in  the  x  direction.  This  increased  the 
resolution  at  the  trailing  edge  and  matches  the  grid  used  in  the 
previous  work  in  which  reversed  flow  was  observed.  This  grid 
was  used  for  the  full  Navier-Stokes  calculation  at  a  >  30  deg 
by  FDL3D-I. 

Each  of  the  simulations  used  the  converged  solution  from 
the  previous  work”  as  an  initial  condition.  The  L2  norm  of 
the  change  in  the  dependent  variables  displayed  a  spike,  which 
then  quickly  decreased  to  its  starting  value.  The  solution  wu 
then  allowed  to  evolve  until  the  L2  norm  reached  and  main¬ 
tained  a  constant  value.  It  should  be  noted  that  the  dependent 
variables  in  ARC3D  are  divided  by  the  Jacobian  whereu  in 
FDL3D-I,  they  are  not.  This  resulu  in  the  L2  norm  differing 
by  several  orders  of  magnitude  between  the  two  codes.  All 
solutions  were  obtained  using  one  processor  on  a  Cray  2 
supercomputer.  The  data  processing  rate  is  expressed  in  ps/- 
point/iteration.  The  data  processing  rates  were  39.46,  43.82, 
and  60.11,  for  ARC3D  laminar,  FDL3D-1  thin  layer,  and 
FDL3D-1  Navier-Stokes,  respectively. 

Resulla  and  Discussion 

The  first  issue  dealt  with  was  that  of  turbulence.  ARC3D, 
with  the  turbulence  model  of  Baldwin  and  Lomax,  wm  run 
I8(X>  iterations  until  convergence.  The  L2  norm  peaked  at 
1 .2  X  Kf*  and  dropped  to  3.9  x  10  ”.  The  average  v^ue  of  y* 
on  the  upper  surface  was  4.01  with  a  maximum  of  9.39  at  the 
trailing  ^ge.  Figure  3  shows  the  pressure  coefficient  distribu¬ 
tion  at  81 . 194k  root  chord  for  the  laminar  and  turbulent  solu- 
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Fig.  3  Cocfflckat  of  prasarc  n  y/T  at  81.11%  foot  chord  for  a 
della  wlag:r>7S  deg;  M.-I.9S;  0-28  deg;  C*-4.4i  x  It*. 


Fig.  4  Coaloan  of  aofaiahwd  pHol  preasare  at  II.lt%  roel  chord 
fora  della  wtag  aalag  ARC3D  wHh  ao  laihoknee  aodel.  T-TS  dag; 
M.-l.»9;a-2tdet,  jr«-4.48  X  It*. 

tions  (u  well  u  data,  which  will  be  discussed  later).  The 
Reynolds  number  based  on  running  length  at  this  streamwise 
location  was  3.63  x  10*.  Also  shown  are  the  data  of  Ouyton’ 
and  that  of  Miller  and  Wood.’  At  this  streamwise  location, 
the  average  value  of  y*  was  3.04,  the  minimum  value  was 
1.76,  and  the  maximum  was  8.70.  The  laminar  solution  fell  be¬ 
tween  the  experimental  resulu  from  the  symmetry  plane  to  the 
secondary  separation  point;  the  turbulent  solution  predicted 
pressure  coeffidenu  that  were  lower  than  the  dau.  The  largest 
difference  between  the  two  numerical  simulations  was  32.8% 
at  the  centerline.  From  the  secondary  separation  point  to  the 
leading  edge,  both  solutions  predicted  lower  pressure  coeffi¬ 
cients  than  the  data.  Near  the  leading  edge,  the  turbulent  solu¬ 
tion  was  12.1%  higher  than  the  laminar  solution,  but  the  dif¬ 
ference  was  less  than  the  scatter  in  the  experimental  data.  It 
appeared  that  the  laminar  solution  compared  more  favorably 
with  the  data. 

Contours  of  the  pitot  pressure  ratio  at  81.19%  root  chord 
are  given  in  Fig.  4  for  the  laminar  solution  and  in  Fig.  3  for  the 
turbulent  solution  The  modeled  turbulence  has  no  effect  on 
the  predicted  location  of  either  the  primary  or  secondsuy 
vortex.  The  effect  of  turbulence  was  only  apparent  near  the 
wall  below  the  primary  vortex  where  the  boundary  layer  was 
thinner  in  the  turbulent  case. 

The  simulated  oil  flow  patterns  for  the  laminar  and  turbu¬ 
lent  solutions  are  presented  in  Fig.  6.  The  primary,  secondary, 
and  tertiary  separation  and  reattachment  lines  are  all  evident. 
They  are,  from  the  leading  edge  inward,  the  primary  separa- 
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Fig.  S  CoBloan  of  ■oranUud  pHol  pnmn  ot  11.11%  root  chord 
for  a  dctta  wlag  aalag  ARCJD  wttk  the  BaMwta-Loaiaa  larbolcocc 
BMdcl:r-7$d(g;M.-1.9S;a>MdcgiJF*-4.4«  x  It*. 


Fig.  t  Slaalalcd  ot  Iracca  for  a  delta  wlag  aalag  ARC3D  tarbakal 
(agger)  aad  laadaar  (lower):  r-7S  deg;  Af. wl.tS;  owM  deg, 
lt*-4M  X  It*. 


tiofl,  lecondary  reaitacliinent,  tertiary  leparation,  tertiary 
reattachment,  iccondary  separation,  and  the  primary  ryat- 
tachment  at  the  centerline.  The  difference  between  (he  laminar 
and  turbulent  shear  stress  patterns  was  very  minor  for  the  en¬ 
tire  surface.  This  indicated  no  major  structural  differences  be¬ 
tween  the  two  simulations.  The  small  changes  in  the  separa¬ 
tion  and  reattachment  lines  were  apparent  in  the  spanwiM 
variation  of  the  shear  stress.  Figure  7  depicts  one  component 
of  the  normalized  surface  shear  streM  r^  u  a  function  of  the 
normalized  span,  at  81.19%  root  chord.  The  shear  stresses 
were  normalized  to  (heir  maximum  at  that  streamwise  loca¬ 
tion.  The  location  of  the  reattachment  and  separation  lines 
corresponds  to  the  change  of  sign  in  t„.  The  secondary  sepa¬ 
ration  wu  slightly  farther  outboard  {(0.27  vs  38.82%  span) 
for  the  turbulent  case,  whereas  the  secondary  reattachment 
wu  inbovd  (89.83  vs  88.84%  span). 

Slight  differencu  in  prusure  coefficient,  pitot  preuure,  and 
shew  strew  distribution  between  the  laminw  and  turbulent 
solutions  were  observed.  The  Baldwin-Lomax  turbulence 
model  is  well  known  to  be  inaccurate  in  regions  of  sepwated 
flow.  For  or  w  20  deg,  where  no  separations  were  seen,  the  ef- 
feas  of  turbulence  were  small.  At  a  -  30  deg,  Iwge  regions  of 
sepwated  flow  have  been  previously  observed,  making  use  of 
(he  Baldwin-Lomax  inappropriate.  Since  we  are  primarily  in- 
teruted  in  the  assewment  of  the  numerical  procedure,  ac¬ 
curacy,  and  flowfleld  structure,  only  laminar  solutions  were 
used  in  the  remainder  of  this  paper. 

The  second  calculation  used  FDL3D-1  with  the  thin-layer 
approximation.  The  solution  required  2400  iterations  using 
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ng.7  StMWWfxsiT- vs  j’/y  at  11.19%  root  chord  for  sMti  wlag: 
r«7S  deg,  Af.>  1.93,0-20  deg,  Jb-4.48  X  M*. 


Rg.  •  Cealean  of  Boraeaitaed  gllol  greteaic  at  11.19%  root  chord 
for  a  deitt  wiag  aeiag  FDUD-l  dda  layer:  r-75  deg,  M.-I.98; 
o-20  4cg;iFe-4.40  x  10*. 

the  tame  time  step  and  damping  coeffleientt  that  were  used  in 
(he  ARC3D  calculation.  The  L2  norm  peaked  at  1 .3  x  (O'’  and 
dropped  to  7.6  X  lO*’.  The  pressure  coefficient  distribution  at 
81.1^  root  chord  is  compared  to  the  ARC3D  thin-layer 
ruult  in  Fig.  3.  The  greatest  difference  between  them  wu 
6.84%  (at  the  centerline)  and,  in  general,  both  were  in  ex¬ 
cellent  agreement  with  the  experimental  data. 

Figure  8  shows  contours  of  the  pitot  preuure  ratio  from 
FDL3D-i  thin  layer  at  81.19%  root  ^ord.  When  compared  to 
the  ARC3D  ruult  in  Fig.  4,  it  wu  seen  that  both  codu  capture 
the  primary  and  secondary  vorticn.  The  predicted  location  of 
(he  primary  vortex  wu  iJmost  exactly  the  same  for  (he  two 
codu.  However,  the  secondary  wu  farther  outboard  for 
ARC3D  (81 .97  vs  84.70%  span).  FDL3D-I  wu  in  better  agree¬ 
ment  with  the  experimental  ruults  from  Monnerie  and  Werle* 
which  placed  the  secondary  at  82.S%  span. 

The  simulated  oil  flow  patterns  compwed  very  well  over  the 
surface  and  are  not  presented  here.  Both  display  the  primary, 
secondary,  and  tertiary  sepwation  and  reattachment  linu.  In 
addition,  both  capture  a  very  small  region  of  reversed  flow  at 
the  trailing  edge,  which  wu  observed  in  the  previous  work. 
The  differencu  between  the  two  thin-layer  solutions  were 
more  appuent  in  the  spanwisc  variation  of  the  shew  strew,  u 
seen  in  Fig.  7.  The  location  of  the  secondary  separation  wu 
newly  the  same  for  both  solutions;  however,  t^  secondary 
reattuhment,  the  tertiary  separation,  and  reattachment  for 
FDL)D-1  were  inbowd  of  their  locations,  u  predicted  by 
ARC3D  (88.02  vs  89.84%  span,  79.62  vs  81.81%  span,  and 


102 


63.38  vs  67.69%  span,  icspeclively).  The  preceding  results  in¬ 
dicate  that  FDL3D-1  (using  the  thin-layer  equations)  and 
ARC3D  produce  very  similar  solutions. 

The  third  calculation  used  the  full  Navier-Stokes  equations 
and  evolved  3000  iterations  with  the  same  time  step  and  damp¬ 
ing  terms  that  were  used  in  the  FDL3D-1  thin-layer  calcula¬ 
tion.  The  L2  norm  peaked  at  l.lSxlO'’  and  dropped  to 
9.7  xIO*’.  The  pressure  coeffident  distribution -at  81.19% 
root  chord  from  the  full  and  thin  layer  Navier-Stokes  equa¬ 
tions  are  compared  in  Fig.  3.  The  two  results  were  nearly  iden¬ 
tical  for  the  first  80%  of  the  span  and  display  a  8.86%  differ¬ 
ence  in  the  pressure  coefficient  near  the  outboard  edge.  There 
was  excellent  agreement  with  the  data  of  Miller  and  Wood.’ 

Figure  9  shows  contours  of  the  pitot  pressure  ratio  at 
81.19%  root  chord  for  the  FDL3D-1  Navier-Stokes  simula¬ 
tion.  Comparison  to  the  thin-layer  result  in  Fig.  8  shows  that 
the  location  of  the  primary  vortex  was  nearly  identical.  How¬ 
ever,  the  secondary  vortex  was  farther  outboard  for  the 
Navier-Stokes  result  (81.97  vs  85.23%  span). 

The  simulated  oil  flows  along  the  upper  surface  again 
showed  only  slight  differences  between  the  two  solutions  and 
are  not  presented  here.  The  shear  stress  from  the  thin-layer 
and  Navier-Stokes  solutions  are  compared  in  Fig.  7.  The 
secondary  separation  point  for  both  results  were  nearly  the 
same,  59.22%  and  58.92%  span  for  the  thin-layer  and  Navier- 
Stokes  solutions,  respectively.  The  secondary  reattachment 


Fig.  9  Cealoars  of  aoraullicd  pitot  pwaw  at  11.19%  root  chord 
for  a  della  wiag  aslag  FOUD-I  Navicr-Slokcs:  r>i7S  dag, 
Af.  - 1 .9$,  a  »  W  deg.  -  4.41  X 10*. 


Fig.  10  Cocffldenl  of  prtasart  vs  y/Y  at  11.19%  root  chord  for  a 
della  wing  aslag  FDUD-I  with  grid  rcfloemcat;  r-7S  deg, 
M.-l.9S;a>20  4eg,  i7«-4.40  x  10*. 


and  tertiary  separation  and  reattachment  all  were  further  in 
board  for  the  thin-layer  result  (92.04  vs  88.02%  spun,  81.72  vs 
79.62%  span,  and  67.87  vs  63.38%  span,  respcciively). 

From  these  above  figures,  it  was  concluded  that  differences 
between  the  thin-layer  and  the  Navier-Stokes  results  were 
quite  small  for  a  delta  wing  at  a  20  deg  and  the  thin-layer 
c^culation  (which  was  less  expensive)  captures  the  essential 
flow  features. 

An  interpolation  using  the  full  Navier-Stokes  solution  was 
then  performed  to  provide  an  initial  condition  for  the  grid  res¬ 
olution  study.  The  fine  grid  had  approximately  twice  as  many 
points  as  the  primary  grid.  In  addition  to  the  doubling  of  the 
computational  domain,  the  local  time  step  was  reduced  by  the 
local  scaling  with  the  Jacobian.  These  two  factors  combined  to 
make  this  an  extremely  expensive  calculation;  thus,  only  600 
iterations  were  run.  This  was  sufFicient  for  the  L2  norm  to 
return  to  its  starting  value.  The  L2  norm  peaked  at  7.03  x  10'* 
and  dropped  to  5.8x10’.  Figure  10  presents  the  pressure 
coefficient  distribution  at  81.19%  root  chord.  There  were  very 
slight  differences  between  the  two  solutions,  particularly  in  the 
region  between  the  primary  and  secondary  vortex.  Contours 
of  the  pitot  pressure  ratio  for  the  fine  grid  solution  are  pre¬ 
sented  in  Fig.  1 1 .  For  the  region  of  the  primary,  secondary, 
and  tertiary  vortices,  the  Tine  grid  solution  reveids  additional 
minute  structures.  These  structures  were  also  indicated  in  the 
lateral  shear  stress  in  Fig.  12.  In  the  fine  grid  soluuon,  t.ie  sep- 


Flg,  II  CoalMFS  of  Bomaliicd  pile*  prtMan  M  11.19%  root  chord 
lor  a  dcMa  wtag  osing  FOUD-I  wttk  grU  rtflocaral:  r-7S  d«g: 
AF..  I.9S;  0-20  drg.  lie -4.41  x  10*. 
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Fig.  12  Sbnr  Hresi  viy/Kti  gl.19%  root  chord  for  a  delta 
aslag  FDUD-I  with  grid  reflaemeal:  F-TS  dag,  AF.-I.95;  o-20 
dag;  ffr-4.4g  x  10*. 
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Fig.  13  Coefflclnt  of  prtMoic  n  jr/V  at  49.M  aad  11.19%  reol 
cfeofd  for  a  Mu  wiog  aalag  FDL3D-I:  T  >>73  deg,  saJS 

deg;  ira-d.4t  X  It*. 
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Fig.  14  Coaloonofa)aoniulliedpHolpmioreaadb)o(/alM.t% 
root  (bord  for  a  ddtt  wlog,  fat  Narler^okca  riaiolatloa;  T  •  7S  deg, 
Af.-l.M;>-)Sdcg;ifa>4.4(  x  It*. 


aration  between  the  secondary  and  tertiary  separation  points 
have  spread  laterally  to  accommodate  the  additional  vortical 
structures. 

The  Navier-Stokes  equations  were  used  to  calculate  the  flow 
over  a  delta  wing  at  a  -  30  deg,  which  required  2000  iterations. 
The  L2  norm  peaked  at  1.2  x  lO*’  and  dropped  to  8.6 x  10-’. 
The  pressure  coefficient  distribution  at  49.80^  and  81.19% 
root  chord  is  presented  in  Fig.  13.  Included  are  experimental 
data  from  Gu^on*  at  SO  and  80%  of  the  root  chord.  From  the 
centerline  to  4S%  span,  the  pressures  were  equal,  but  out¬ 
board  they  diverge  with  the  downstream  pressures  being  up  to 
S.0%  higher.  The  agreement  with  the  experimental  data  wax 
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Fig.  W  SlaaMcd  od  Iraca  for  a  defta  wlag  adag  Ihte  layrr  tapper) 
aad  Navtcr-Slokes  Oower);  r>79  deg;  M^^l.95;  a-33  deg: 
ife-4.a  X  10*. 


very  good.  Figures  14  show  contours  of  normalized  pitot 
pressure  (Fig.  14a)  and  the  product  of  the  density  and  the 
streamwise  Mach  number  (Fig.  14b)  at  89.74%  root  chord. 
This  product  was  negative  in  the  darkened  region,  indicating 
reverted  flow.  The  reverted  flow  wu  contained  in  a  region 
along  the  surface  beneath  the  secondary  vortex.  Results  from 
the  previous  work  are  presented  in  Figs.  IS.  The  location  of 
the  primary  vortex  it  the  same  for  the  two  simulations.  How¬ 
ever,  the  h^ght  of  the  reverted  flow  region  (normalized  by  the 
span)  was  much  smaller  (0.0884  compared  to  0. 19S).  TTiis  dif¬ 
ference  is  due  to  the  additional  diffusion  calculated  in  the 
Navier-Stokes  solution. 

Figure  16  depicts  the  surface  shear  map  for  the  present 
Navier-Stokes  solution  and  the  previous  thin-layer  solution. 
The  reversed  flow  region  extends  forward  to  46%  root  chord 
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compared  to  for  the  thin-layer  solution.  Also  visible  are 
two  large  counter-rotating  vortical  structures.  In  the 
nomenclature  of  Tobalc  and  Peake'*  these  are  referred  to  as 
nodes.  These  vortices  are  formed  by  the  viscous  interaction  be¬ 
tween  the  reversed  flow  and  the  streamwise  flow.  The  inboard 
vortex  was  swept  up  and  inward  by  the  primary  vortex.  The 
outboard  vortex  was  swept  up  and  inward  by  the  secondary 
vortex.  These  two  struaures  merge  above  the  reversed  flow  re¬ 
gion  and  are  swept  downstream.  In  the  Navier-Stokes  solu¬ 
tion,  the  inboard  and  outboard  nodes  are  located  at  84. 13  and 
8S.I^  of  the  root  chord,  respectively.  In  the  thin-layer  solu¬ 
tion,  they  are  located  at  M.23  and  76.92%  of  the  root  chord, 
respectively.  In  addition  to  the  large  nodes,  a  smaller  node  was 
located  in  each  of  the  solutions.  These  nodes  are  too  small  to 
be  discerned  in  Fig.  16.  For  the  Navier-Stokes  solution,  the 
third  node  is  located  at  53. 1  %  root  chord  along  the  secondary 
separation  line.  In  the  thin-layer  solution,  the  third  node  is 
located  at  the  trailing  edge  along  the  tertiary  separation  line. 

These  are  significant  differences  between  the  thin-layer  and 
Navier-Stokes  simulations  for  a  delta  wing  at  a  =  30  deg. 
These  differences  manifested  themselves  in  the  secondary 
vortex  where  a  shear  layer  is  present.  The  variation  in  the 
streamwise  and  normal  (to  the  upper  surface)  extent  of  the 
reversed  flow  region  is  most  evident. 

Condusloiu 

The  supersonic  flowfleld  about  a  delta  wing  (P=75  deg) 
was  numerically  simulated  at  M.  =  1.95  and  4.48  x  10*.  Two 
different  codes  based  on  the  approximately  factorized  implicit 
Beam-Warming  algorithm  were  used  in  the  study.  To  assess 
the  effecis  of  turbulence,  ARC3D  was  used  with  and  without 
the  Baldwin-Lomax  turbulence  model.  The  laminar  and  tur¬ 
bulent  solution  exhibited  very  little  difference  in  the  predicted 
vortical  structure.  This  indicates  that  turbulence  was  not 
overly  critical  for  a  delta  wing  at  a  =  20  deg  with  the  present 
grid  resolution. 

The  comparison  between  ARC3D  and  FDL3D-1  (thin  layer) 
showed  minor  differences  in  the  predicted  pressure  coefficient 
distribution,  location  of  the  primary  and  secondary  vortices, 
and  surface  shear  stress.  A  comparison  was  made  b^ween  the 
thin-layer  and  Navier-Stokes  solutions  generated  by  FDL3D-I. 
The  two  solutions  were  nearly  identical,  indicating  that  there 
was  no  advantage  in  using  the  Navier-Stokes  equations  over 
the  thin  layer  equations  of  delta  wing  simulations  at  moderate 
angles  of  attack. 

A  calculation  carried  out  on  a  refined  grid  for  a  =  20  deg 
reveals  that  the  numerical  result  is  essentially  independent 
from  grid  refinement.  A  calculation  with  FDL3D-I  using  the 
full  Navier-Stokes  equations  was  then  made.  For  the  a  •  20  deg 
calculation,  there  was  no  improvement  with  grid  resolution. 

Navier-Stokes  solution  for  a  delta  wing  at  a =30  deg 
showed  an  embedded  reverse  flow  region  that  was  observed  in 
an  earlier  thin-layer  calculation.  The  dimensions  of  this  sepa¬ 


rated  flow  region  was  significantly  different  from  the  thin- 
layer  calculation.  The  reversed  flow  region  exhibited  strong 
velocity  gradients  in  directions  other  than  normal  to  the  body 
surface.  The  thin-layer  approximation  does  not  account  for 
the  viscous  term  associate  with  these  gradients. 
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